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Abstract. High order and fractional PDEs have become prominent in theory 
and in modeling many phenomena. Here, we focus on the regularizing effect 
of a large class of memoryful high-order or time-fractional PDEs — through 
their fundamental solutions — on stochastic integral equations (SIEs) driven by 
space-time white noise. Surprisingly, we show that maximum spatial regular- 
ity is achieved in the fourth-order-bi-Laplacian case; and any further increase 
in the spaiiaZ-Laplacian order is entirely translated into additional temporal 
rcgularization of the SIE. We started this program in [T][5|, where we intro- 
duced two different stochastic versions of the fourth order memoryful PDE 
associated with the Brownian-time Brownian motion (BTBM): (1) the BTBM 
SIE and (2) the BTBM SPDE, both driven by space-time white noise. Under 
wide conditions, we showed the existence of random field locally-Holder solu- 
tions to the BTBM SIE with striking and unprecedented time-space Holder 
exponents, in spatial dimensions d = 1,2,3. In particular, we proved that the 
spatial regularity of such solutions is nearly locally Lipschitz in d = 1,2. This 
gave, for the first time, an example of a space-time white noise driven equation 
whose solutions are smoother than the corresponding Brownian sheet in either 
time or space. 

In this paper, we introduce the 2/3 _1 -order /3-inverse-stable- Levy-time Brow- 
nian motion (/3-ISLTBM) SIEs, p e {l/2 fc ;fceN}, driven by space-time 
white noise. Based on the dramatic regularizing effect of the BTBM den- 
sity (/3 = 1/2), and since the kernels in these /3-ISLTBM SIEs are fundamental 
solutions to higher order Laplacian PDEs; one may suspect that we get even 
more dramatic spatial regularity than the BTBM SIE case. We show, how- 
ever, that the BTBM SIE spatial regularity and its random field third spatial 
dimension limit are maximal among all /3-ISLTBM SIEs — no matter how high 
we take the order 1/p of the Laplacian. This gives a limit as to how far we 
can push the SIEs spatial regularity when driven by the rough white noise. 
Furthermore, we show that increasing the order of the Laplacian P~ 1 beyond 
the BTBM bi-Laplacian manifests entirely as increased temporal regularity of 
our random field solutions that asymptotically approaches that of the Brown- 
ian sheet as p ~\ 0. Our solutions are both direct and lattice limit solutions. 
We treat many stochastic fractional PDEs and their corresponding higher or- 
der SPDEs, including BTBM and /3-inverse-stable- Levy-time Brownian motion 
SPDEs, in separate articles. 
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1. Introduction, motivation, and statement of results 

Lately, many phenomena in mathematical physics, fluids dynamics and turbu- 
lence models, mathematical finance, and the modern theory of stochastic processes 
have been related to and described through deterministic fractional and higher order 
evolution equations (e.g., see 014], [6]-[n], [19], [23]-[26], [31], [34]-[37], 001 E], 
[43]-[48], and [53]); and it is only natural to investigate these important equations 
under the influence of a driving random noise. 

In the two recent articles [T] [5] we introduced two new stochastic versions of 
fourth order memory-preserving (which we coin memoryful) deterministic PDEs 
related to Brownian-time processes (BTPsjQ — introduced in [Mj [13] — driven by 
space-time white noise: 



A BTP, in its simplest form, is a process X x (\Bt\) in which X x is a Markov process starting 
at x S K d and B is an independent one dimensional BM starting at 0. A Brownian-time Brownian 
motion (BTBM) is a BTP in which X x is also a Brownian motion. BTPs include many new and 
quite interesting processes (see 14 13 31 147|). which we are currently investigating in several 
directions (e.g., [31 M Q3])- With the exception of the Markov snake of Le Gall ([42]), BTPs 
fall outside the classical theory of Markov, Gaussian, or semimartingale processes. We label BTP 
PDEs as memoryful since the initial data is part of the PDE itself (see JL4}) 
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(1) the space-time- white- noise-driven Brownian-timc Brownian motion (BTBM) 
SPDE 

(d t U=^ + ±A 2 U + a(U)dl+ d W, (i,i)e(0,»)xR d ; 
[U(0,x) =u (x), xeM d , 

where dl~^. d W is the space-time white noise on R + x K d — and on a proba- 
bility space (fi, J^",P) — that corresponds to the Brownian sheet W; and 

(2) the stochastic integral equation we called BTBM SIE 

■'BTBM" „ . /„ \ i . . / / tt^BTEM 1 * 



(1.2) U(t,x)= K»™My)dy+ K?™ y a(U(s,y)W(dsxdy) 

JR d JR d JO 

where K^.™ 1 is the density of a <i-dimensional Brownian-time Brownian 
motion given by: 

/>oo 

(1-3) K*f= 2 K^ y K™ s ds 

Jo 

with K™". u = e ~,'*~wa 2 ' and = ^-J^; and where is the white 

noise on R + x R d . 

Unlike the deterministic case a = 0, (jl.ll) and (|1.2[) behave differently, and each 
is quite interesting in its own right. Each of these two equations gives a different 
stochastic interpretation of the memoryful BTBM PDE in [14j [13] : 

( i4) l dtU = 7Bi + l A2u; {t > x) e (0 ' TO) x K "' 

[u(0,x) = uo(z); xel d . 
and its equivalent integral for nfl 

(1.6) u(t,x) = [ K*™ d u (y)dy. 

jR d 

As proven in [51 [T], the SIE (|1.2p — which we also denote by e|™ BM (a,-u ) — nas 
real random field solutions in d = 1,2,3 with striking Holder regularity in which 
the time-space Holder exponents are ^^jp , (pp A l) J , as we recall precisely 



2 For a review of the BTPs higher order and fractional PDEs connections and generalizations, 
as well as connection to the important Kuramoto-Sivashinsky PDE, we refer the reader to |14l 
US] [6) US] H3 148] and the references therein. The connection of BTPs to their fourth order PDEs 
(including JT3}) was first given in |14| . Also, their connection to time-fractional PDEs was first 
established implicitly via the half derivative generator in 1 141 . In 46 47, 48 the equivalence 
between a large class of high order and time-fractional PDEs, including il .41 1 and 

t 1 1 

d?u= — Au; t £ (0,oo), a; S M d , 
(1.5) { 1 VS 



u(0,x) = uo(x); 

was established explicitly, using the Caputo fractional derivative. For a discussion of interesting 
aspects of these PDEs see also the introduction in [T| . In the new multiparameter-time case the 
reader is referred to l4ll3l. 
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in Section 11.11 belown, and it is similar in regularity to the following L-Kuramoto- 
Sivashinsky (L-KS0SPDE 

(17) id t U = -±A 2 U-±AU-±U + a(U)dt+ d W, (t,x) G (0, oo) x R d ; 

1 U(0,x) = u (x), lel^, 

obtained from the linearized KS PDE in [6] by adding a multiplicative space-time 
white noise term (see [2]). In [9j [10], we treat a large class of higher order and 
fractional — and rougher — SPDEs, including (jl.ll) and its equivalent time-fractional 
SPDE 

{dfu=-j=AU + a(U)dl+ d W; t G (0, oo), x G R d , 
\t/(0,a;) = u (x); 

i 

where <9 t 2 is a factional derivative in time (see e.g. |43] ) . 

In this article, we focus on a large class of fascinating stochastic integral equa- 
tions driven by space-time white noise and generalizing the BTBM SIE (|1.2[) : the 
/3-inverse-stable-Levy-time Brownian motion SIEs (/3-ISLTBM SIEs), which we dis- 
cuss in more details in Section [fT2"l below. These SIEs are obtained from the BTBM 
SIE in (11.21) by replacing the BTBM density with the fundamental solution to the 
2/3" 1 = 2u order, /3 _1 € {2 fe ; k € N}, memoryful PDEs 

(1 . 9) -g^J^ + ^S^; M^MX* 

(up(0,x) =u (x), x G R d 

and their equivalent time-fractional PDEs 

(1 1Q s I flf «/»(*, a:) = ^A^(t, a:) (t, a;) G (0, oo) x R d 

I tt^(0,x) = uo(af), a; G M d , 

whercQ Ep tK — S^&I^ ; the process is the /3-inverse-stable-Levy motion de- 
scribed in Section [TT21 below, and 5f is the well known Caputo fractional derivative 
of order (3 G {l/2 fc ; k £ N} in time (see e.g. [43]). 

Based on the dramatic regularizing effect of the BTBM density on the space-time 
white noise driven BTBM SIE (|1.2[) as just described above (see also Theorem ll.il 
below), and due to the fact that the kernels in the /3-ISLTBM SIEs of this article are 
fundamental solutions to the higher order PDEs (|1.9p ; one may suspect that we get 
even more dramatic spatial regularity than the BTBM SIE case, possibly obtaining 
random field solutions in arbitrarily high spatial dimensions as j3 \ (y /*• oo) 
instead of just d = 3 as in the BTBM case (/3 = 1/v = 1/2). We show, however, 
that the BTBM SIE spatial regularity and its random field third spatial dimension 
limit are maximal among all /3-ISLTBM SIEs; no matter how small we take /3 (how 



^In particular, as was established in [2], the BTBM SIE 11.21 1 has nearly locally Lipschitz 
solutions in d = 1,2. This fact provided for the first time a counterexample to the common 
folklore non-wisdom that "a solution to a space-timc-white-noise-driven equation cannot have a 
solution that is more regular, temporally or spatially, than the Brownian-sheet in the underlying 
white noise" 

4 The L in the name refers to the linearized PDE part. Such L-KS SPDE is treated in [2]. 
^As usual, E denotes the expectation operator. 
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high we take the order /3 _1 of the Laplacian). Further, we show that increasing the 
order /3 _1 of the spatial Laplacian beyond the BTBM order of 2 translates entirely 
into temporal regularization of our /3-ISLTBM SIEfl This surprising result is the 
regularity content of our two main theorems: Theorem [L2] and Theorem 11.31 below . 

1.1. Recalling the Brownian-time Brownian motion SIE case. Before stat- 
ing our first main result, it is instructive to recall the BTBM SIE results in [Tpl . 
Following pQ, we denote b}@ fF*~'^(T x R d ;R) the space of real-valued locally 
Holder functions on Txl d whose time and space Holder exponents are in (0,7i) 
and (0,7 S ), respectively. The first main result in [1] is now restated. 



Theorem 1.1 (Allouba 1 ). Fix < 7 < 1. Assume the following Lipschitz 
and growth conditions 

{(a) \a(u)-a(v)\ < C \u - v\ ii,»eM; 
(b) a 2 {u) <C{l + u 2 ); 
(c) M € C^ 7 (R d ;R) and nonrandom , V d = 1,2,3. 

hold. Then there exists a pathwise-unique strong solution (U, W) to 
e S B T BM (a,Uo) on R + x R d , for d = 1,2,3, which is L P (Q.) -bounded on T x R d 

for all p > 2. Furthermore, U G H 8 A 2 J (T xR d ;R) for every 
d= 1,2,3. 



Theorem 11.11 states that the stochastic kernel integral equation (|1.2|) has ul- 

3 — 

tra regular strong/] solutions on R+ x R d , namely U € Hs (T x R; R), U € 

1- 1- 1- 

H4 (T x R 2 ;R), and U € Hs '2 (TxR 3 ;R). I.e., in space, the BTBM paths 

have a rather remarkable — and initially-surprising — nearly local Lipschitz regular- 
ity for d = 1,2; and nearly local Holder 1/2 regularity in d = 3. This is remarkable 
because the BTBM kernel is able, in d = 1, 2, to spatially regularize such solutions 
beyond the traditional H61der-l/2~ spatial regularity of the underlying Brownian 



"i.e., the extra regularizing "energy" of spatial Laplacians of orders higher than that of the bi- 
Laplacian is converted to extra temporal regularity, when faced with the extremely rough driving 
space-time white noise. 

7 Earlier, in [5], the additive noise case a = 1 for e g™ BM ( a > u o) was considered; and the existence 
of a pathwise unique continuous BTBM SIE solution U(t, x) for x £ R d and d = 1,2,3, such that 



sup E P |l/(t,:r)| 2p < C 



(4-d)p 
l+t * 



t > 0, d = 1, 2, 3, p > 1, 



was proved. 

8 Throughout the paper, T = [0, T] for some fixed but arbitrary T > 0. Here and in the 
sequel C^' 7 (K <i ;R) C Cj(K d ;R) denotes the space of bounded p-times continuously diffcrentiable 
functions such that all derivatives up to (and including) the p-th order are bounded and all p- 
th order derivatives are Holder continuous, with some Holder exponent < 7 < 1. Also, the 
boundedness conditions on uq and its derivatives may easily be relaxed as in [3]. 

^Here strong is in the stochastic sense of the noise if and its probability space (f2, & ', {J^},P) 
being fixed a priori. Throughout this article, whenever needed, we will assume that our filtrations 
satisfy the usual conditions without explicitly stating so. 
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sheet corresponding to the driving space-time white noisd^l This degree of smooth- 
ness is unprecedented for space-time white noise driven kernel equations or their 
corresponding SPDEs; and the BTBM SIE is thus the first such example. In time, 
our solutions are locally 7-Holder continuous with dimension-dependent exponent 
7 G (0, ^jp) for d = 1,2,3. This is in sharp contrast to traditional second order 
reaction-diffusion (RD) and other heat-operator-based SPDEs driven by space-time 
white noise, whose fundamental kernel is the Brownian motion density and whose 
real- valued mild solutions are confined to the case d — 1. In this regard, the di- 
chotomy between the rougher paths of BTBMs as compared to standard Brownian 
motions on the one hand (quartic vs. quadratic variations) and the stronger reg- 
ularizing properties of the BTBM density vs. the BM one on the other hand is 
certainly another interesting point to makJ^l. 

1.2. The /3-inverse-stable-Levy-time Brownian motion SIE: the first main 
theorem. In the first main result of this article, we generalize the first BTBM SIE 
result in [1] Theorem 11.11 to the interesting case of the inverse-stable-Levy-time 
Brownian' motion SIE with index /3 = l/u, v e {2 fc ;fceN} (/3-ISLTBM SIE), 
which we now motivate and introduce^- This generalization allows us to better 
appreciate how hard it is to smooth away space-time white noise. 

1.2.1. Recalling (3-ISLTBM. Inverse stable subordinator — which we also call 0- 
inverse-stable-Levy motion and denote by — arise in the work of Meerschaert 
et al. [45j 04] as scaling limits of continuous time random walks. Let S(n) — 
Yi + • • • + Y n a sum of independent and identically distributed random variables 
with EY n = and EY% < 00. The scaling limit C -^ 2 S{[ct}) B(t) as c — > 00 
is a Brownian motion B at time t, which is normal with mean zero and variance 
proportional to t. Consider Y n to be the random jumps of a particle. If we impose 
a random waiting time T n before the nth jump Y n , then the position of the par- 
ticle at time T n = J\ + • • • + J n is given by S(n). The number of jumps by time 
t > is N(t) — max{n : T„ < t}, so the position of the particle at time t > 
is S(N(t)), a subordinated process. If P(J„ > t) = t~&l(t) for some < /3 < 1, 
where l(t) is slowly varying, then the scaling limit c _1 /^Tj cf ] =>■ Lp(t) is a strictly 
increasing stable Levy motion Lp at time t and with index (3, sometimes called a 
stable subordinator. The jump times T n and the number of jumps N(t) are inverses 
{N(t) > x} = {T(\x~\) < t} where \x] is the smallest integer greater than or equal 

"As noted in pQ, it is important to note here that the common "folklore wisdom" of solu- 
tions of space-time-white-noise driven equations not being smoother than the associated Brow- 
nian sheet — in either space or time — originated from the predominant case of SPDEs, in which 
either the underlying kernel is that of a Brownian motion or the spatial operator is a Laplacian. 
The kernel ^ , however, is much more regularizing to the space-time-white-noise driven 
e BTBM ( a ' n °) than the density of BM is to its corresponding equation. This becomes evidently 
clear in Lemma l2.4l Lemma l2.3l and Lemma l2. 21 (compare to the more traditional BM and random 
walk case in |12p. 

^We observe in passing that — roughly speaking — the paths of e S? BM ( a i u o) in d = 1 are 
effectively 3/2 times as smooth as the RD SPDE paths in d = 1, in d = 1 the BTBM SIE is as 
smooth as an RD SPDE in d = 1, and in d = 3 our BTBM SIE is half as smooth as an RD SPDE 
in d = 1. Also, for d = 2, 3, the spatial regularity is roughly four times the temporal one, and in 
d = 1 the spatial regularity is maximized at a near Lipschitz vs near Holder 3/8 in time. 

12 Throughout this article we assume that u = f}~ 1 £ {2 fe ; k £ N}, where N is the set of natural 
numbers. The case /3~ 1 = 2 is the BTBM SIE case, with a minor scaling of the Brownian motion 
as discussed in [3]. 
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to x. It follows that the scaling limits are also inverses c~@N(cb) A^(t) where 
Ap(t) = ini{x : L(x) > t}, so that {Ap(t) < x} — {Lp(x) > t}. We call the process 
a /3-inverse-stable-Levy motion. Since N(ct) « c^A(i), the particle location may, 
for large c, be approximated by C -^/ 2 S(N([ct])) « (c' 3 )- 1 / 2 5(c^A (3 (t)) w B(Ap(t)), 
a Brownian motion subordinated to the inverse or hitting time (or first passage 
time) process of the stable subordinator hp. The random variable Lp{t) has a 
smooth density. For properly scaled waiting times, the density of Lp(t) has Laplace 
transform e~ tsP for any t > 0, and the random variables Lp(t) and t^-^LpiX) 
are identically distributed. Writing gp(u) for the density of Lp{l), it follows 
that Lp(t) has density t~ x I ^ g p{tr x I ^ u) for any £ > 0. Using the inverse rela- 
tion F(Ap(t) < x) — f{Lp{x) > t) and taking derivatives, it follows that Ap(t) has 
density 

(1.11) K&„ = t^x-^gpitx-^), 

As noted above, we assume throughout this article that = /3 _1 G {2 ; fe G N}. 
In this case, there is a simple connections between fc-iterated Brownian-time Brow- 
nian motion and /3-ISLTBM. We denote by 

M\ (t):=B*(\B k (..-B 2 (\B 1 (t)\) ■■■)]) 

OS. 

a fc-iterated Brownian-time Brownian motion at time t; where {Bi}^ =1 are inde- 
pendent copies of a one dimensional scaled Brownian motion starting at zero, with 

density -7==^ exp (~|t ) = (l/^)^™,^ anc ^ independent from the standard 

d-dimensional Brownian motion B x , which starts at 2: e R d . By B^ fc (t) = 

B x (A 1 / 2 fe(^)) we mean a d-dimensional /3-ISLTBM — with /3 = l/2 fe — starting at 
lei 11 and evaluated at time t; in which the outer BM B x and the inner A 1 / 2 k are 
independent. 

Lemma 1.1 (The /3-ISLTBM density). The probability distributions ofW k (t) 

OB. 

i=l 

and B^ k (t) are the same for every k = 1,2,... and every t > 0. In particular, 
when /3 — l/2 k , k £ N, the Ap and the (3-ISLTBM transition densities are given by 



(1.12) k -2 



respectively 



Proof. Let /3 = l/2 fe , fc e N. By Corollary 3.1 in g7] we get that the distributeions 

1 We are using the convention n<=!o c ' = 1 f° r an y c i ana - * ne convention Jjg f(s)ds = /(s), 
for every /. Also, we use the convention that the case k = (/? = 1) in the /J-ISLTBM is the 
standard d-dimensional Brownian motion case. 
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are the same. Now, equation (0.14) in [3] gives us that 

(1.13) ^ = ^ exp (_f!)=_L^. 

This, together with Lemma 3.1 and Lemma 3.2 in |47j and a simple conditioning 
argument using the independence of all the Brownian motions, we immediately 
obtain (|1.12[> as asserted. □ 

We now define our /3-ISLTBM SIE as the stochastic integral equation: 

(1.14) Uf,(t,x)= [ K™^u (y)dy + [ f ^?.£>a(U p {a, y))W{ds x dy) 

JR d JR d JO 

where K^ 1 y Af> is the transition density of a <i-dimensional /3-ISLTBM, starting from 
x G R d , B x Aft := {B x {Kp{t)),t > 0} given bj0 

(1-15) K^/' = / K™ y Kilsds. 

Jo 

We also denote the /3-ISLTBM SIE (PI) by e^f gLTBM (a, u Q ). Just as in the BTBM 
SIE case, e ^ I f SLTBM ( a J u o) is one of t w0 different stochastic versions^ of the higher 
order (2v = 2/3 -1 ) memoryful PDEs (jl.9p and their equivalent time fractional PDEs 

Of course, in the deterministic case, both (|1.9p and (|1.10l) are equivalent to their 
integral form 

(1.16) u p (t,x)= [ K^dy. 

Jm d 

1.2.2. First theorem: 2/3 _1 order SIEs regularity and third dimension maximality. 
Our first main theorem is now stated. 



Theorem 1.2 (Spatio-temporal regularity and third dimension maximality: 
direct solution). Fix f3 = X/v, v £ |2 fc ;fc £ N}. Assume the following Lips- 
chitz, growth, and initial smoothness conditions 

'(a) \a(u)-a(v)\ < C \u - v\ u,v £M.: 
(Lip) I (b) a 2 (u) < C(l +u 2 ); u £ K, 

(c) u a G C^ _2 ' 7 (R d ;R) and nonrandom , V d= 1,2,3. 

/jo/gL TTjen i/iere exists a pathwise-unique strong solution {U$,W) to 
e SIE (a, u ) on M.+ x R d , for d=l,2. 3, which is LP(U)-bounded onTxR d 

' %V-d\ ~ /" 4-d ^ }/ 



/or aZ/p > 2. Furthermore, Up £ H V 4,y / ' V 2 ' (Tx M d ; K) /or even/ 
d= 1,2,3. 



14 Compare with the expression of ' A<3 Lemma 11.11 in terms of scaled BM transition 



densities 
15 Thi 

11.9) or from l|1.10|l by adding the white noise term as in [9] 



^The other stochastic version is the 2i> or the time-fractional f3 order SPDE obtained from 
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Theorem O states that, for /3 = 1/u and v € {2 fc ;fceN}, these 2/3" 1 or- 
der /3-ISLTBM SIEs have quite interesting locally-Holder solutions with temporal 
and spatial Holder exponents given by ( 2 ^~° ) and A l) , respectively, for 
d = 1, 2, 3. Comparing this regularity with the corresponding result for the fourth 
order BTBM SIE in Theorem ll.il we see that the spatial regularity — spatial Holder 
exponent and the maximum spatial dimension of 3 — is identical. Since, the fun- 
damental density (fundamental solution) estimates leading to the regularity con- 
clusions of Theorem 11.21 — Lemma 12.21 to Lemma 12.41 — are sharrPl this means that 
there is a limit as to how far we can push against the powerful roughening effect of 
the driving space-time white noise. Despite the fact that these SIEs are co-driven by 
fundamental solutions of arbitrarily high order (2f3~ 1 ) PDEs involving the spatial 
/3 _1 -Laplacian operators, we can obtain locally Holder real random field solutions 
only up to three spatial dimensions and with spatial Holder exponents up to the 
maximal BTBM bi-Laplacian case (/3 _1 = 2), for all v = (3~ x e {2 fe ;fc e N}, no 
matter how large /3 _1 is. 

To appreciate the richness of the regularizing effect of these /3-ISLTBM SIEs, 
however, we need to look beyond just the spatial dimensionality and regularity 
aspects. So, we will now examine the conclusion of Theorem 11.21 regarding the 
maximum temporal (effective) Holder exponent, as /3 \ 0. As observed above, the 
strong roughening influence of the space-time white noise prevents further spatial 
smoothing of our /3-ISLTBM SIEs beyond the BTBM bi-Laplacian case, no matter 
how large /3 _1 gets. However, all of the extra smoothing "energy" resulting from 
increasing the spatial Laplacian order /3 _1 cannot simply be "destroyed" by the 
white noise; and it is converted instead into temporal regularization of these /3- 
ISLTBM SIEs (as /3 \ 0). Theorem 11.21 describes precisely this temporal effect in 
terms of Holder exponents. In particular, the maximum effective regularity of the (3- 
ISLTBM SIEs increases asymptotically to the well-known Holder (l/2)~ regularity 

of the Brownian sheet; i.e., the maximum effective Holder exponent ^ 2| ^_7 d ) /* 

\ as (3 \, for every d = 1,2,3. The following table summarizes our regularity 
findings and compares them to the more standard and classical case of reaction- 
diffusion SPDEs driven by space-time white noise. 



d 


Random Field Solutions 


Holder Exponent (time, space) 




RD SPDE 


/3-ISLTBM SIE 


RD SPDE 


/3-ISLTBM SIE 


1 


Yes 


Yes 


(i ) 




2 


No 


Yes 


N/A 




3 


No 


Yes 


N/A 


((^r,(§r) 



Table 1.1. /3-ISLTBM SIEs (y = /3" 1 e {2 fc ;fceN}) vs RD 
SPDEs (/3=1). 



16 We will have more to say about the regularity of these /3-ISLTBM SIEs in We also 

briefly note that by third dimension maximality, we mean maximality among integer dimensions. 

^The effective Holder exponent is the minimum of the spatial and temporal Holder exponents, 
which of course determine how smooth the random field solutions arc as functions of both time 
and space together. 
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To prepare for the statement of our results under the less-than-Lipschitz condi- 
tions in QNLipD (Th eorem [TT3l below) . we now introduce the spatial lattice version 
of e SIE (a, tin) as well as introduce the new associated process we call /3-inverse- 

/3-ISLTBM v ' u/ 

stable-Levy-time random walk and define the lattice limit solutions involved in the 
statement of Theorem 11.31 The main machinery we use in the proof in this case is 
our K-martingale approach, which we introduced and used in the BTBM SIE case 
in [I]. We recall this approach, adapting it to our settineFj in Section [4721 

1.3. The spatial lattice version and the second main result. As in pQ, we 

now spatially discretize e ^ I f SLTBM ( a ; u o)- This accomplishes at least two things: 
(1) it gives a multiscale view of e ^ : f SLTBM ( a i u o) an d (2) it allows us to prove our 
existence and regularity results without the Lipschitz condition on a. 

1.3.1. (3-inverse-stable-Levy-time random walk on the lattice. In [181 112] . standard 
continuous-time random walks on a sequence of refining spatial lattices 

J X£ :=[]{■■•> ~ 26 ^ °> S n, 2«n, ■ ■ •} = 5 n lA 

I i=l ) n>l 

(with the step size S„ \ as n oo) played a crucial role — through their 
densities — in obtaining our results for second order RD SPDEs. In [T], in the 
fourth order Brownian-time setting, that role is played by Brownian-time random 
walks on X^: 

(1.17) S x B<5n (t):=Sl(\B t \); 0<t<w,xeX* 

where Sg (t) is a standard rf-dimensional continuous-time symmetric RW starting 
from x £ X^ and B is an independent one-dimensional BM starting at 0. The 
subscript 5 n in (|1.17|) is to remind us that the lattice step size is S n in each of the 
d directions. 

In this article, we replace Brownian-time random walk with /3-inverse-stable- 
Levy-time random walk (/3-ISLTRW): 

(1.18) S% e<Sn (t)~Sl(A (t)); 0<f3< 1, 0<t < oo, x e X^ 

It is then clear that the transition probability (density) 5 ™ ' ff of the /3-ISLTRW 
S^ itfB (i) onX^ is given bj0 

d a r°° d 

(1.19) K^ 5 ;' fi = 2 / K^yK^ s ds; < < 1, < t < oo, x,y eXf 

Jo 

where K t . x ™ is the continuous-time random walk transition density starting at 
x £ X^ and going to y € X^ in time t, in which the times between transitions are 

exponentially distributed with mean S^f. I.e., K^ x Sn is the fundamental solution 
to the deterministic heat equation on the lattice X^ : 

(1.20) _2i2 = _a„<(;) ; (f,s)e(0,oo)xX* 



18 A11 we need to adapt it here is to replace the BTRW kernel of [T] with the /3-ISLTRW one 

in tl. 191 below. 

in KWjf ,Ab KWf ,As 

Throughout this article, K. t ™ := K. t . x " (with a similar convention for all transition 

densities). 
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where s^ n '■= A„/2 is the generator of the RW Sg n (t) on X^. 

By mimicking our proof of Theorem 0.3 in 0] , we easily get a 2v order differential- 
difference equation connection to /3-ISLTRW: 

Lemma 1.2 (/3-ISLTRW's DDE). Fix [3 = l/v, v e {2 k ;k e N}. Let u^ n (t) = 
E 



u ^ j (i)J with uo as in ( |NLip[ ) . Then up^ n solves the following 2v order 
differential- difference equation (DDE) on R+ x X^ : 

du^t) = ^ A^) 1 

(1.21) I dt ^ 2 K t 1 ~ K / u p ' 2^ ^' 
< l 

,«|,„(0) = uo(a:)i 

w/iere E PiK = Ef(A ^ (1)) " . Moreover, K^*"'^ so/ues (TOTl) on [0,oo) x X^ ; with 

(1.22) u a (x) = K . x = K Q . X = 

1.3.2. Lattice (3-ISLTRW SIEs and limits solutions to /3-ISLTBM SIEs. The crucial 
role of the /3-ISLTRW density in our approach to the /3-ISLTBM SIEsCEH]) becomes 
even clearer from the following definition of our approximating spatially-discretized 
equations: 

Definition 1.1 (Lattice /3-ISLTRW SIEs). By the /3-ISLTRW SIEs associated 
with the BTBM SIE e™ M (a,Mo) we mean the system {e^ SLTRW (a, u , n)} 
of spatially-discretized stochastic integral equations on R + x X^ given by 

\ - ^Rwf .A, _ , , , \ - f*Twi,As.~ v ,^dWy(s) 




(i.23) uut) = j2 KS'^ofe) + e / 



where the /3-ISLTRW density is given by (I1.19[) . For each n g N, we think of 
{W*(t);t > 0} as a sequence of independent standard Brownian motions indexed 
by the set X^ (independence within the same lattice). We also assume that if 
m ^ n and x G Xf„ n X^ then W£{t) = W%(t), and if n > m and x e X^ \ X^ then 
W£(0 = 0. 

Notation 1.1. We will denote the deterministic and the random parts of (|1.23j) by 

Uf <nD {t) and Up^ n ^{t) (or U^ D {t) and Ug° R {t) when we suppress the dependence 
on n), respectively, whenever convenient. 

We define two types of solutions to /3-ISLTRW SIEs: direct solutions and limit 
solutions. 

Definition 1.2 (Direct /3-ISLTRW SIE Solutions). A direct solution to the f3- 
ISLTRW SIE system {e^ SLTRW (a,u ,n)j on R + x with respect to the 
Brownian (in t) system {W£(t);t > 0}^ n I ) gMxX d on the filtered probability space 
(O, J^", {^t},P) is a sequence of real-valued processes < U n > with continuous 

I J n— 1 

sample paths in £ for each fixed x € X^ and n € N such that, for every (n, x) € 
N x Xfj, Up n {t) is ^-adapted, and equation (|1.23[) holds P-a.s. A solution is said 
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to be strong if {W£(t); t > 0} („ )X )gNxX <i ancl (f2j IP) are fixed a priori; and 

with 

(1.24) J^ = o-{<r(W^(s);0<s<t,a;eX£,neN)u<yr}; t £ R+, 

where ./f is the collection of null sets 

{O : 3 G G Sf , O C G and P(G) = 0} 

and where 

yt>o 

A solution is termed weak if we are free to choose (O, {^t}, IP) and the Brownian 
system on it and without requiring & t to satisfy (|1.24[) . Replacing M + with T := 
[0, T] — for some T > in the above, we get the definition of a solution to the 

/3-ISLTRW SIE system <^ e SIE (a, u , n) } on T x R d . 

f J |_ ,3-ISLTRW v ' u ' y J n=1 

The next type of /3-ISLTRW SIE solutions we define is the first step in our 
K-martingale approach of [I] , which we recall in Section 14.21 By first reducing 
e ^ I f SLTRW ( a j u o> Ti) to the simpler finite dimensional noise setting, it takes full ad- 
vantage of the notion of /3-ISLTRW SIEs limit solutions to /3-ISLTBM SIEs. 

Definition 1.3 (Limit /3-ISLTRW SIE Solutions). Let I £ N. By the ^-truncated 
/3-ISLTRW SIE on K+ x X^ we mean the /3-ISLTRW SIE obtained from {T23} by 
restricting the sum in the stochastic term to the finite <i-dimensional lattice X„ ; := 
X^ n {[— I, l] d ; I £ N} and leaving unchanged the deterministic term Up n D (t): 

t 



(1.25) U'l^t) 
where 



ui n . D (t) + E / ^ s ,t(^ (s) ) dw "™ (s);xeX " 



u 



(Ulrjr)) := -J^a(Uy^(r)), Vr,s < t. 

On 

We denote (|1.25|> by e *^ TRW -(a» w °' n ' ^ x n ^ ^' a sonl ti on t° tne system of 
truncated /3-ISLTRW SIEs I e*"^ TKW .(o, u , ra, J) [ on R+ x X^ with respect to 
the Brownian (in t) system {W%(t);t > 0} xeXd on the filtered probability space 
(f2, {^ti, IP) is a sequence of real-valued processes \ Up n l \ with continuous 

sample paths in t for each fixed x £ X^ and I £ N, such that, for every (l,x) £ 
N x Xf i; Up nA {t) is ^-adapted, and equation (fl~25|) holds P-a.s. We call t/^„ 
a limit solution to the /3-ISLTRW SIE fl~2"3]) if Up, n is a limit of the truncated 
solutions Up^ n j (as I — > oo). If desired, we may indicate the limit type (a.s., in L p , 
weak, . . . , etc). 



So, by 



Remark 1.1. In both (Ti~25l) and (Ti~23l) . Ul n D {t) = E u Q feaj*)) 

Lemma [1.21 Up nD (t) is differentiable in time i and satisfies (|1.21j) . Also, using 
linear interpolation, we can extend the definition of an already continuous-in-time 
process Up n {t) on K + x X d , so as to obtain a continuous process on K + x R d , for 
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each n G N, which we will also denote by Up n (t). Henceforth, any such sequence 
{Up^n} of interpolated Up^'s will be called a continuous or an interpolated solution 
to the system I eS pf SLTRW ( a i u 0: n ) r • Similar comments apply to solutions of the 
truncated e t_SIE (a.Un.n.l). 

/3-ISLTRW v ' u ' ' ; 

We now define solutions to e ^ I f SLTBM (o, uq) based entirely on their approximat- 
ing < e SIE (a.un.n) >, through their limit. Since we defined direct and limit 

b l /3-ISLTRW v ' u ' ' J ' b 

solutions to e ^ I f SLTRW ( a i u o, n), for each fixed n, we get two types of /3-ISLTRW 
SIEs limit solutions to e ^f SLTBM ( a ! u o)- direct /3-ISLTRW SIEs limit solutions and 
,5-ISLTRW SIE double limit solutions. The "double" in the second type of solutions 
reminds us that we are taking two limits, one from truncated to nontruncated fixed 
lattice (as I — > oo) and the other limit is taken as the lattice mesh size shrinks to 
zero (as S n \ or equivalently asn/ oo). 

Definition 1.4 (,3-ISLTRW SIEs limits solutions to e s Jf arrr ,(a, u )). We say that 
the random field U is a /3-ISLTRW SIE limit solution to e SIE (a, tt ) on K+ x R d 

r /3-ISLTBM v ' + 

iff there is a solution {£/? „(i)}neN to the lattice SIE system \ L (a, u , n) \ 

on a probability space (f2, {^t},P) and with respect to a Brownian system 
{W%(t);t > 0} (n^jgNxx 1 sucri that J7 is the limit or a modification of the limit 

of \Ug n f (or a subsequence thereof). A /3-ISLTRW SIE limit solution U is 

called a direct /3-ISLTRW SIEs limit solution or a /3-ISLTRW SIEs double limit 
solution according as {Ug n (t)} n ^ is a sequence of direct or limit solutions to 

{ e ^ I fsLTRw^ a ' UQl n )} ' r ^^ LC um ^s may be taken in the a.s., probability, L p , 
or weak sens^3- We say that uniqueness in law holds if whenever and 
are /3-ISLTRW SIEs limit solutions they have the same law. We say that path- 
wise uniqueness holds for /3-ISLTRW SIEs limit solutions if whenever {^n^} an d 

are lattice SIEs solutions on the same probability space and with respect 
to the same Brownian system, their limits and are indistinguishable. 

1.3.3. Second main theorem: the lattice-limits solutions case. We can now state 
our second main result of the paper. The following theorem gives our lattice- limits 
solutions result for e ^ I f SLTBM ( fl j Uq) under the non-Lipschitz conditions ( |NLip[ ) on a. 

Our limits solutions result under Lipschitz conditions is stated in Theorem I A . lIFI . 
which is proved in Appendix [AJ 



When desired, the types of the solution and the limit are explicitly stated (e.g., we say 
strong (weak) /3-ISLTRW SIEs weak, in probability, L P (Q), or a.s. limit solution to indicate that 
the solution to the approximating SIEs system is strong (weak) and that the limit of the SIEs is 
in the weak, in the probability, in the L P (Q), or in the a.s. sense, respectively). Of course, we may 
also take limits in any other suitable sense. 

21 The type of limit solutions in the Lipschitz case is direct limit solutions as opposed to the 
double limit solution in Theorem 1 1.3 1 
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Theorem 1.3 (Spatio-temporal regularity and third dimension maximality: 
lattice-limits solutions). Fix ft = \ jv. v € {2 fe ;/c G N}. Assume the conditions 

(a) a(u) is continuous in u; u£l. 
(NLip) I (6) a 2 {u) < (7(1 + u 2 ); u E R, 

(c) uq G C^ I/_2 ' 7 (K d ; R) and nonrandom , V d = 1,2,3. 

hold. Then, there exists a /3-ISLTRW SIE double weak-limit solution to 
ef® SLTBM {a,u ), U, such that U(t,x) is L p (0,P) -bounded on T x R d for ev- 

' 2v — d\ ~ ( i—d ^ i 



en/ p>2 and Up e H V / V 1 ) (fx R d ; R) /or every d = 1, 2, 3 



Remark 1.2. Of course, we can use change of measure — as we did in our earlier 
work on Allen-Cahn SPDEs and other second order SPDEs (see e.g. jTTl [THl US] and 
all our change of measure references in [T^] for results and conditions) — to transfer 
existence, uniqueness, and law equivalence results between e ^ : f SLTBM ( a , uq) and the 
/3-ISLTBM SIE with measurable drift e SIE (a, b, u Q ): 



(1.26) 



Up (t, x) = I u (y)dy + I f K™*.^ KUp (s, y))dsdy 

R d JO 



H-t£>(U P (s,y)W(ds x dy), 



under the same conditions on the drift /diffusion ratio. If it is desired to inves- 
tigate e ?f SLTBM ( a i b, uq) on a bounded domain in M. d with a regular boundary, 

we simply replace the /3-ISLTBM density K™ J/ ' A ' i! in (j 1 . 26|) with its boundary- 
reflected or boundary-absorbed version (the /3-ISLTBM density in which the outside 
d-dimensional BM is either reflected or absorbed at the boundary). 

The proof of Theorem 11.31 under the conditions ( |NLip[ ) is neither standard nor 
straightforward — even after obtaining the new non-trivial spatio-temporal regular- 
ity estimates (in Lemma 12.31 and Lemma 12.41 below) on the unconventional kernel 

I^t;x y' Af3 • This is because standard techniques, like the classical martingale prob- 
lem approach, do not apply directly to kernel equations like the /3-ISLTBM SIE 
e^_f SLTBM (a,Mo) or its discretized version e^ SLTRW (a, u , n) under flNLipD . This 
leads us to use our K- martingale approach, introduced in [1]. 

2. Key estimates 

2.1. Density regularity estimates and third dimension maximality. The 

first set of estimate^ we need are bounds on the square of the /3-inverse-stable- 
Levy-time Brownian motion density ^ A/3 and its associated lattice /3-inverse- 

stable-Levy-time random walk density K^ 5 "' ? and their temporal and spatial 
differences. We obtain these estimates for both kernels simultaneously. The method 



22 As is customary, all constants may change their value from one line to the next without 
changing their notation. Also, to simplify notation, we will often suppress the dependence on /3 
without further notice. We will denote the Euclidean norm on d-dimensional spaces by |-|. 
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of proof is to reduce, via an asymptotic argument, these estimates for the /3-ISLTRW 
to the corresponding ones for the /3-ISLTBM density ' 3 and perform the 
computations in the continuous setting of the /3-ISLTBM. Since all the results in 
this part hold for all n > N* (equivalently for all 5 n < 5n*) for some positive integer 
N* , we will suppress the dependence on n, except when it is needed or helpful, to 
simplify the notation. Also, whenever we need these estimates, we assume that 
n > N* without explicitly stating it every time; and when we do, we letF^l 



(2.1) 



:= {n EN;n> N*} 



We start by observing that in the classical setting of Brownian motion and 
its discretized version continuous-time random walk on X^ = 5 n 2i d , we have the 
following well known asymptotic result relating their densities (see e.g., |52| ) 



(2.2) 



Kf™* y S d as n oo (as S n -> 0); Vi > 0, x,y £ 



^HSn >[»]«»! 

where for each x £ M. d we use [x]s n to denote the element of X^ obtained by 
replacing each coordinate x^ with S n times the integer part of S~ Xi, and a n ~ b n 
as n — > oo means a n /b n — > 1 as n — > oo. Now, for every continuous and bounded 
d — > K, we have 



(2.3) lim £ 

3/GX^\{x} 



«o(y)^ 



r ™/^o(y)<%; t > o, x e M d , d > l, 



and by the dominated convergence theorem we obtain 



lim 

S„\0 



(2.4) 



E 



*;Ws„,Mf. 



wo(y) 



E 

yGX^\{:c} 



uo(i/)<# 



lim > 

2/eX£\{x} 



if" 



for t > 0, x € R d , and d > 1; since, by (|2~2]l 



lim ^ _ftT ' , 1" r , 

yex^ 



«o(y) 



lim 

<5„V) 



My) > ^t ; o )S rfs 



,.u (y)dy 



for every (s,x) € (0, oo) x M. d . We then straightforwardly get the following result. 



Lemma 2.1. For every continuous and bounded uq 



and for every d > 1 



(2.5) lim x 

yex^ 



^ X f K^ A ^ (y)dy;V(t, a; )e(0,oo)xR rf , 



and the following asymptotic relation holds between the (3-ISLTBM and (3-ISLTRW 
densities: 



(2.6) K 



r BM d ,Ap <. d 
H;x,y °n 



as n — y oo (as S n — > 0); t > 0, x,y £ R d , x =^ y. 



2 ^We adopt these simplifications with lattice computations throughout the paper. 
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Remark 2.1. Equation ([2. 5ft confirms the intuitively clear fact that the kernel 
form of the /3-ISLTRW DDE (|1.21j) converges pointwise — as S n \ — to the kernel 
form of its continuous version, the /3-ISLTRW PDE in |3] . We also remind the 
reader that the right hand side of ((23j) is in C 1 ' 2 ^ 1 for all (t, x) £ (0, oo) x R d 
under the uq conditions in ( |NLip[ ). 

Our first regularity lemma for the densities is now stated. It implies, among 

other things, that there is a considerable smoothing effect of K™ ' Af> as /3 gets 
smaller; however it also implies that our SIEs don't possess random field solutions 
beyond the third spatial dimension, no matter how small /3 gets. 

Lemma 2.2 (Smoothing and third dimension maximality). There are constants C 
and C , depending only on d and /3 = 1/v, v £ |2 fe ; k £ N}, and a S* > such that 
for all S < 8* 



r BM a ,Ap 
H;x 



dx = Ct 2u and 



2 _ -d 

< C5 d t2l7; 



for all t > 0, d — 1, 2, 3. Hence, 







fl 




/0 JR d 





' 2v~d 

dxds = Ct 2v and 



xer" 



! . , 2v-d 

ds < C8 d t~^r-. 



for allt > 0, d = 1, 2, 3. In addition, J R 
oo, /or aZZ d > 4. 



r BM ,A^ 
4;x 



Jo /k 



rBM ,Kp 
^s-x 



dxds 



Proof. First, fix an arbitrary f3 1 = v £ |2 fe ,fceN}. Using the definition of 
Kj.^ ' A|9 , Lemma [2.11 and Lemma [1.11 here together with Lemma 3.1 and Lemma 
3.2 in S3 we obtairJ3 



(2.7) 



lim > 

x£X d 



<5 d 



rBM ,Ay3 



OO /"OO 



Jo 

OO />OG 



8i;x U\;x UjU ' 



JO 

oo />oo 



JO 



1 

[2tt(si +ui)] 
2 k 



d/2 



Kt-is^fAu.dsidm 



K t . 

(O.oo)*- 1 ' v3 



[27r(« 1 +ui)] ,,/a 

fe-2 



fc-2 

I'.M TT r/BM 

afc-<;0,- 

i=0 



_! o!s2 • • • (is ^ 



V5 



21 



Recall that we are using the convention J* r0 f(s)ds = f(s), for every /. 
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Gathering the two inside integrals and transforming to polar coordinates (s,, Uj) i— > 
(pi, #i), letting p — (pi, . . . ,pk) and = (6±, ... ,6k), and noticing that all pi for 
i = 2, 3, . . . ,pk cancel when k > 2; equation (|2.7p become^ 



C 



2 fe-2 r pi_,-i^ 2 (f fc -i-i) , pi_i-i"" 2 Cfc-»-i) 



fc-2 



(2.8) (o, 7 r/2) fc (o,«,)» p| ^ [sin(^) + cos(0 x )]* Vsin(0fc- l )cos(0 fc _ l ) 



i=0 



1(7*2*,; d=l,2,3, 

|oo; d > 4. 
Then there is a S* > such that, whenever 5 < 6*. we obtain 



J. A, 



< Ct2v - : d= 1,2,3, 



with a finite constant C > C. The last assertion of the lemma trivially follows upon 
integration over the time interval (0,t]. □ 

The following lemma is key to our Holder regularity result in time. We give a 
probabilistically-flavored proof using the notion of 2-/3-inverse-stable-Levy-timcs 
random walk and 2-/3-inverse-stable-Levy-times Brownian motion given below. 

Lemma 2.3 (Kernel temporal regularity). There is a constant C , depending only 
on d and (3 = v G |2 fc ;fc G N}, and a S* > such that for 5 < 5* 



(2.9) 




S~— s;x 



2 2i/-d 

dxds < C(f - r) 21/ , 

2 , 2i/-d 

ds < CS d (t-r)^T, 



for < r < t and d = 1,2,3, with the convention that K™*' Af3 = = K|™ M<i £/ 
t < 0. 



2 ds<C5 d {t-ry^r- d = 1,2,3. 



Proof. We will prove that 

(2.io) r £ 

for all £<<$*, for some 5* > 0, simultaneously with its corresponding /3-inverse- 
stable-Levy-time Brownian motion density statement. The first step is to show the 
identity 

2 



(2.11) 



RW?,A 



(t— r)\x 



*a+(t—r),a+(t—r);0 ' 



2K 



RW| n ,2As 
s+(t-r),s;0 



25 Equation l2~8t is the reason for the third spatial dimension maximality. 
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where 



(2-12) 



K 



u.v.Q 



OO /"00 



HVV 5 „A0 J . 

JX r 1 +r 2 ;0 JX u;0,r 1 rL v;0,r2 U ' lUI 2 



is the density of the 2-/3-inverse-stable-Levy-times random walk 



(2.13) § A « A c*>,JM : = (A^(«) +Af («)) ; < u,v 



< oo, 



in which the d-dimensional random walk S® (on X^) and the two identically- 
distributed one-dimensional processes and A^ are all independent. But, 



\ " lir HW|,A^ llf HW? ) Ap 



(2.14) 



OO pOO 



^0 



E 



K tfi, ri K ts),r 2 dr idr2 



oo /.oo 



JQ 



Jx r 1 +r 2 ;0 ri -u;0,r 1 r> -v,0,r 2 u ' l a/ : 



r RW d n .2A/3 



2 - ^u,v-0 



The identity (|2.11l) immediately follows from (I2.14j) . Similarly, we get the corre- 
sponding identity for the /3-inverse-stable-Levy-time Brownian motion setting 



(2.15) 



dx 



rBM _ 1K -BM d ,A^ 

\+{t-r);x ^ s ' x 

,BM d ,2A 3 K -BM d ,2A p _ glK -BM*,2A^ 

S+(i-r),s+(t-r);0 + ^sjO /Jt S+(i-r),s;0 



where 



(2.16) 



rBM ,2A^ _ 



OO /'OO 



JO 



K BM a K ^ii TC fin rfro 

A ri+r 2 :l) n ii;0,ri n »;0J2 ( "l , "2 



is the density of the 2-/3-inverse-stable-Levy-times Brownian motion 
(2.17) X° A(1) A{2) (u,v) :=X°(A«(n) + Af («)) ; 0<n, V <oo, 



in which the d-dimensional BM A" and the two identically-distributed one-dimensional 
processes A^ and are all independent. Using the identities (|2.11l) and (|2.15l) . 
along with a similar asymptotic argument to the one we used in the proof of 
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Lemma 12.21 together with the dominated convergence theorem, yield 



lim • 



(2.18) 



t 



r RW a n ,2Ap 
^s+(t-r),s+(t-r);0 



ds 



= lim / y 

J0 xex* 



S+(t-r);x ^ s > x 



-ds 



ds-2 



^s+(t-r),s:0 



ds 



Kit 



A(i _ w BM d ,Ap 
{t-r);x ^ s ' x 



dxds 



BM«,2A^ s _ 2 



4 d r* 

K+{l-r),a+{t-ry,O ds + / K ~s~s;6 
"0 

t ft ft 

]K 2s+ 2(t-r)ds + / K 2s ds - 2 / K 2s+ ( 4 _ r )ds 
o Jo Jo 



,BM d ,2A^ 



ds 



t-r 
2 



^ 2s ^s - / K 2s ds - 

I t — T 



/+ 



lisds 



2t-r 



ds 



t+ 



for d = 1,2, 3, where K w is defined in terms of K^,.Q 2A/i by the relation 



(2.19) K -- K «.-;o 

(it, u) = (a, a) or (it, v) = (a + 6, a) or (it, i?) = (a, a + 6) for some a, b > 0. 



u; = u + u and (it, v) has one of the forms 



We observe that 



K 2u = lC ;0 2A ^ = / / K^ r2:0 KXX-kr 2 dndr 2 



OO /"OO 



JO 

oo c OO 



(2.20) 



o Jo 



K BM K BM dx 



rBM ,A« 



da; = CM2l7; d= 1,2,3 



The last assertion follows from the computation in (|2.7p and (|2.8[> . It is clear then 
that K 2u is decreasing in u, for every v = 1/(3 g {2 fc ;fce N}. Thus, the sum of the 
last three terms of the (|2~T8)) is < 0. This and ([2^20)) give us (|2~T0|) for all S < S* , 
for some 5* > and for some constant C > 0, together with its corresponding 
/3-inverse-stable-Levy-time Brownian motion density statement; and Lemma 
follows at once. □ 



The following spatial difference second moment inequality for the /3-ISLTRW and (3- 
ISLTBM densities reflects their critical spatial-regularizing effect on our solutions. 
The following lemma captures the surprising fact that we cannot improve on the 
spatial regularity of the BTBM SIE by decreasing f3 below 1/2. This implies the 
maximality of the BTBM SIEs spatial regularity among the family of /3-ISLTBM 
SIE family. 
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Lemma 2.4 (Kernel spatial regularity). Let (3 G {l/2 k ;k & N} and define the 
intervals 

C (0,1]; d=l, 
h =1(0,1); d = 2, 
l(0,|); d = 3. 

for any given positive numbers {ad G -^d}d = i; t^ere exists a constant C depending 
only on (3, d and {ad}^ =1 , o,nd a 5* > such that for S < 5* 



(2.21) 



JR d L 

/SI 



RWf ,A 3 _ ^RWf 



"dxds < C\z\ 2ad t pi - ad ^\ 
2 ds < C5 d \z\ 2ad t p{ad ^\ 



for t > 0, where < C < oo and < p(ad, /?) < 1 /or every ad G /or d = 1,2,3 
and /or every /3 G {l/2 fc ; k G N}. 

Remark 2.2. For a given /3 _1 G {2,3,4, . . .}, and on any compact time interval 
T = [0, T], the inequality (12.21)) may — for any given value ad — be rewritten as 



(2.22) 




r BM ,Ap 



~ds < CS d \z\ 2ad ; 



dxds < C\z\ 2a "; 



where, for each d = 1,2,3 



C = C sup T^ a "-^ < oo 

a d <^I d , 

/3e{i/2 fc ; /ceN}. 



also depends on T in (12.22[) . 



Proof. Let j3 = 1/2 for k G N. Starting with the L 2 estimate involving the spatial 
difference of the /3-ISLTBM density in (|2.21j) . letting u\ — using the polar 
transformation (r^, ui) i-> (pi, letting p = (pi, . . . , p^) and £ = . . . , 0f-), and 
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(2.23) 



rBM ,Aj tti-BM 



dxds 



r-t pOO pOO p J' 

/ / / / U( K ?£ -^MX^dndr, 
Jo Jo Jo Jwt d Z.1 v ' 



ds 



t coo coo 



10 JO JO 

rt coo coo 



wBM J _ojf BM<i W AjS K Af3 dridind* 
ZJX r 1 +u 1 :0 ZJX r 1 +u 1 ;z I A s;0,n A s;0,«i" r l™l ttS 



1 _ g 2(r 1 +u 1 ) 

/o 7o ./o [27r(ri + ui)] d ^ 2 

fc-2 \ 



(0,oo)* 



'75 r <=- 
i=0 

fc-2 



n/2 



/ ^Co 2| TT o a=i=i d "2 ■••***) driduids 

7(0,00)*-! ' £_„ u *-<- '— 73— / 



< c 



2 k ~ 2 _ f fk-i-1 

I - r ] e ^ e L 4pfc - 

i=0 



fe-2 



-dpdOds 



< C 



(0,f )* (o,«.)» p= X s [ sm (0) + cos(0)] * J] A /sin(0fc-i)cos(0 fc _ i ) 

i=0 

a fe-2 fl-i-i l 



1 - e 2 fi e 



i=0 



(0,oo) fc 



I- 1 

Pi s 



-dpds 



z\ 2a e— 



< C 



f 






oo) k 


Pi 


2a^p 1 (a.j3) 


, d = 


2a-f-p 2 {a,l3) 


d = 




d = 



fe-2 

II' 



-dpds 



for some finite constants Cj, z = 1,2,3, where C2 and C3 depend on <o, and 
where we have used the simple facts that min <e<7r/2 [sin(0) + cos(0)] = 1 and that 
1 — e~ u < u a for u > and < a < 1. This proves the L 2 estimate for the 
/3-ISLTBM density in (|2.21l) . Then, an asymptotic argument similar to the one in 



26 See Remark [2721 in [T] for a detailed discussion in the BTBM case ,8=1/2. 
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the proofs of Lemma 12.21 and Lemma 12.31 yields 
(2.24) 




together with the desired /3-ISLTRW density L 2 estimate in (|2~2T|) for all 5 < S* , 
for some S* > 0, with possibly different constants. □ 



2.2. Spatio-temporal estimates for /3-ISLTRW and /3-ISLTBM SIEs. In 

this subsection, and assuming only the less-than-Lipschitz conditions ( |NLip[ ) on 
a — together with a temporary moment condition — we obtain spatial and temporal 
differences moments estimates that are crucial in obtaining the regularity of the 
/3-ISLTRW SIE e^_f SLTRW (a,w ,n) for each fixed n G N* (see (|2~T]) L the tightness 

of the /3-ISLTRW SIEs sequence \ e^ IE (a,u ,n) \ , as well as the Holder 

I /3-ISLTRW x J n£N' 

regularity for their limiting /3-ISLTBM SIE. To make it more convenient for the 
proof of our first main result in the direct solution case, Theorem II .2i we include 
the corresponding spatio-temporal statements for the /3-ISLTBM SIE in the same 
lemmas, together with those for their lattice cousins. 

Fix n G N*, and assume Up n solves e™ 5 (a, u , n) in (|1.23j) and Up solves 
the /3-ISLTBM SIE e^f (a, u Q ) in (|1.14|) . Suppressing the dependence on n, 
let Mp, q (t) = sup x E\U^{t)\ 2q , and M 0>q {t) = sup K E\Up(t, x)\ 2q for q > 1 and 
/3 € {l/2 fe , k e N}. Writing Up and Up in terms of their deterministic and random 
parts Up{t) — Up D (t) + Up R (t) and Up(t, x) = Up t o{t, %) + Up t B,{t, x), we observe 
that Up D (t) is smooth in time by Lemma ll.2l and Up_o(t,x) is smooth in time and 
space as it is a solution to PDEs of order 2/3 -1 as in [31 07] • The next two lemmas 
give us estimates on the random part. 

Lemma 2.5 (Spatial differences). Assume that ( |NLip[ ) holds and that Mp^ q {t) and 
Mp, q {t) are bounded on any time interval T = [0,T]. Th ere exists a constant 
C depending only on q > 1, max x | wo (a:) 1 7 (3 = v G |2 fe ;fc € N}, the spatial 
dimension d — 1, 2, 3, ad, and T such that 

2q 



(2.25) 



\2qa d 



U x p, R (t) - U y p R {t) <C\x~y 
Up, R (t, x) - Up, R (t, y)\ 2q < C\x - y\ 2 ^ 



for all i,j/£ K d , t € T, and d = 1,2, 3; where {a>d} d=1 are as in Lemma \2Ji\ I.e., 
in d = 1, we may take otx = 1; in d = 2 we may take any fixed ct2 € (0, 1); and in 
d — 3, a3 may be taken to be any fixed value in (0, 4). 



Proof. We prove the lattice SIE statement in (|2 .25[) for f/^g; the proof of the 



2 This is the aforementioned temporary moment condition. It is assumed here (in Lemma l2.5l 
and Lemma |2 . 61 below) only to simplify the presentation and to get to the proof of Theorem 1 1.1 1 as 
quickly as possible in Section [3] In Section |4. II this moment condition is shown to automatically 
hold under l |NLip[ l. 
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statement for Up follows the exact same steps, with obvious modifications and will 
not be repeated. Using Burkholder inequality, we have for any (t, x, y) G T x X 2d 



(2.26) 



■2q 



< CI 



E 

zGX d 







H-s:x,z ^t-s;y,z 



ds 



For any fixed but arbitrary point (t, x, y) G T x X 2d let ^ ' v be the measure defined 
on [0, t] x X d by 



d/J-t (s,z) 



rZ.s 



and let \(J%' y \ = /^ :W ([0,i] x X d ). We see from (|2.26[) . Jensen's inequality applied 
to the probability measure Ht' v / \p, t |, the growth condition on a, the definition of 
Mp tq (t), and elementary inequalities, that we have 



uin® - ui R (t) 



2? 



< CE 

< c 



[0,t]xX d 



a(U$(s)) 



i x,y\Q 



1 + Mp,g( S ) 



(2.27) 

< c\ 

l[0,t]xX d 

Now, using the boundedness assumption on Mp_ q on T for d = 1, 2,3, we get 
E 



2q 



Up*® ~ ^ c k e T ^ [^ Pd(Qd) J I* - y\ 2qad 

<C d \x-y\ 2 « a *; a d el d , 

where the last inequality follows from Lemma 12.41 and (|2.22[) in Remark 12.21 and 
where the constant C < oo is as in Remark 12.21 □ 



Lemma 2.6 (Temporal differences). Assume that ( |NLip[ ) holds and that Mp t q(t) 
and M/3 t q(t) are bounded on any time interval T = [0, T]. There exists a constant 
C depending only on q > 1, max x |uo(x)|, (3 = v E |2 fe ;fc G N}, the spatial 
dimension d = 1, 2, 3, and T such that 



(2.28) 



E\Up, R (t,x)-Up >R {r,x)\ 
E Ul R {t)-Ul R {r) ^ <C\t 



{2v-d)q 

2q <C\t~r\~^r- 



r\ 

{2v-d)q 



x G R d ,t,r G T, 
x G X d ,t,r G T, 



for d= 1,2,3. 



Proof. We prove the lattice SIE statement in (|2.28[) for Up; the proof of the 
statement for Up follows the exact same steps, with obvious modifications. Assume 
without loss of generality that r < t. Using Burkholder inequality, and using the 
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change of variable p = t — s, we have for (r, t, x) € T 2 x X rf 
(2.29) 



2<; 



< CE 



■CI 



E 
E 



t— s;x,z r — s;x 7 z 



2 r7s 



For a fixed point (r, £, i) and a fixed /3, let p,% t r be the measure defined on [0, r] x X d 

by 

j X ( \ r„,HwlA« „RwlA fl l 2 rfs 



H— s;x,z r— s;x,z 



5 d 



and let tr \ = fJ-p t r ([0, r] x X d ). Also, for a fixed a; e X d and /3, let k x be the 
measure defined on [0, t - r] x X d by 

2 dp 



<frsf (p) = 



6 d 



and let = K%([0,t — r] x X d ). Then, arguing as in Lemma l2~5l above we get that 

E 



Uln(t) - Ul R (r) ** < C (| M |, t , r | 9 + < C(< - rf^ , 



for <i = 1,2,3, where the last inequality follows from Lemma [2.21 and Lemma [273 
completing the proof. □ 



3. Proof of the first main Theorem 

Here, we prove Theorem 11.11 We start first by recalling a useful elementary 
Gronwall-type lemma whose proof can be found in Walsh |54j . 

Lemma 3.1. Let {<?n(i)}^L ^ e a sequence of positive functions such that go is 
bounded on T = [0,T] and 

9n(t)<C f g n ^(s)(t~s) a ds, 71=1,2,... 
Jo 

for some constants C > and a > — 1. Then, there exists a (possibly different) 
constant C > and an integer k > 1 such that for each n > 1 and t G T 

/"* f-s 

9n+rnk(t)<C m g n (s)- — G?S; 771=1,2,.... 

Jo (m-1)! 
We are now ready for our proof. 

Proof of Theorem ] For the existence proof, we construct a solution iteratively. 
So, given a space-time white noise W, on some (fi, J^", {J^" t }, P), define 



vf\t,x) = / K^u (y)dy 

t 

U^ X \t,x)=Uf\t,x) + f f KTS;xf y a(U { n \s,y))W(d S xdy) 



m d Jo 
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We will show that, for any p > 2 and all d = 1, 2, 3, the sequence I Ug^ (t, x) > 

I p J n>l 

converges in L p (fl) to a solution. Let 



D 0>n>p (t,x) := E Uj" +1) (t,x) - Uffcx) 



D p,n,p(t) SU P Dfs^^t^x). 

Starting with the case p > 2, we bound Dp <niP using Burkholder inequality, the 
Lipschitz condition (a) in ( |Lip[ ), and then Holder inequality with < e < 1 and 
Q =p/(p-2) to get 



< C 

< C 



E 



t—s;x,y 



a{U^\s,y)) - a(U^ L '( S ,y))\ W(ds x dy) 
U^{s,y)-U^{s,y)' 



An-l), 



m d Jo 
t r 



dsdy 



R d JO 



R d Jo 
t 



r BM ,Ap 
^t—s;x,y 



2eq \P /2q 

dsdy I 



[ K t-s-x,y) Df}, n _i tP (s,y)dsdy 



Take e = (p — 2)/p in the above (2eq = (1 — e)p = 2), take the supremum over the 
space variables, and use Lemma |2~21 to see that, for d = 1,2,3 the above reduces to 



p-2 



(3.2) 



^,n,p(*) < C U — 



The case p = 2 is simpler. We apply Burkholder's inequality to D n: 2 and then take 
the space supremum to get 



(3.3) 



D}, n , 2 (t)<C D; >n _ lt2 (s)[t-s}2u ds 



I.e., on any time interval T = [0, T], the integral multiplier on the r.h.s. of (|3.2j) is 
bounded; and if Dp n _ x is bounded on T then so is D% , for every p > 2. Now, 



D} fi p (t) < C sup E 







I f 


¥ BM i ,A (J ' 






jR d Jo 





(U { °\s,y))dsdy 



Since uq is bounded and deterministic, then so are U^ 1 and a(U^). The latter as- 
sertion follows from the growth condition on a in ( |Lip[ ). Thus, by Lemma l2.2l _D^ 
is bounded on T for d = 1, 2, 3 and so are all the „ p . Lemma f3.ll now implies 

i/p 



that for each d = 1,2,3, the series X)m=o 



(t) converges uniformly 

i/p 



on compacts for each n, which in turn implies that Y^=o n converges 
uniformly on compacts. Thus Up converges in L p (fl) for p > 2, uniformly on 
T x R d for d = 1,2,3. Let Up(t,x) :— lim n _).oo ui n \t,x). It is easy to see that [/^ 
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satisfies (flTTljl . and hence solves the ,3-ISLTBM SIE e^ SLTBM (a,u )< This follows 
from (|3.1[) since the Lipschitz condition in (|Lip[) gives 



o(C^(t s a;))-a(C^ n) (t,a:)) 



< CI 



Uf,(t,x)-Ujp(t,x) 



— > as n — > oo 



uniformly on T x l 11 . Therefore, the stochastic integral term in (|3.ip converges to 
the same term with f/g replaced with the limiting Up — i.e., it converges to the 
corresponding term in e^f (a, uo) — as n 

Id - I i-i X J_5 1V1 



E 

< C 



K d JO 



oo, for 

f-i^ (a(C^(*,»)) ~ a(^ n) ( S ,t/))) ^(cfe x 



-i 2 



/ /' 


'BM d ,Af3~ 


2 

E 


m d Jo 







Up{s,y) -Up(s,y) dsdy 



as n 



oo. It follows that Up satisfies the /3-ISLTBM SIE e^ 1 ® f a , uq). Also, 

the solution is strong since the Up are constructed for a given white noise W, 
and the limit Up satisfies (|1.2[) with respect to that same #\ Clearly Up is L p (57) 
bounded on T x R d , d = 1, 2, 3, for any p > 2 and for any T > 0. 

To show uniqueness fix an arbitrary /3 _1 € {2 fc ; fc € N} — and suppress the de- 
pendence of solutions on j3 — and let d = 1, 2, 3, let T > be fixed but arbitrary, and 
let U\ and U 2 be two solutions to our /3-ISLTBM SIE (|TT4|) that are L 2 (0)-bounded 
on T x R d . Fix an arbitrary (t,x) e M+ x R d . Let D(t,x) = U 2 (t,x) - Ui(t,x), 
L 2 (t,x) = ED 2 (t,i), and L* 2 (t) — swp xeRd L 2 (t, x) (which is bounded on T by 
hypothesis). Then, using (jl.141) . the Lipschitz condition in ( |Lip[ ), and taking the 
supremum over the space variable and using Lemma 12.21 we have 

rt 



L 2 (t,x) 



M d Jo 



HU 2 (s,y)) - aiU^y))] 



r BM d ,Ap 
H— s\x,y 



dsdy 



(3.4) 



< C 



R d Jo 
t 



L 2 {s,y) 



<C L*(s) 



K 



^t—s;x,y 

, - i 
BM ,A^ 

t— s:x.y 



dsdy 



dyds < C 



(t- S)2u 



-ds 



Iterating and interchanging the order of integration we get 




L* 2 (r) 



L2(s)ds 



ds 



_d_ d_ 

(t - s) 2v (s - r) 2v 



dr 



L* 2 (t)<C[ / V 2 {s)ds 



for every t > 0. An easy application of Gronwall's lemma gives that L* 2 = 0. So for 



every (t, x) € 



l d and d — 1, 2, 3 we have U\(t, x) — U 2 (t, x) with probability 



one. The indistinguishability of U\ from U 2 , and hence pathwise uniqueness, follows 
immediately from their Holder regularity, which we now turn to. 
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For any given 8 1 = v 6 \2 k ;k 6 N}, we have just shown that, under the 
Lipschitz conditions QLipD , our /3-ISLTBM SIE in ifO^l has an L p (ft)-bounded 
solution Up(t,x) on T x M. d for any T > and any p > 2. Equivalently, Mp iQ (t) = 
snp x E\Up(t, x)\ 2q > 1, is bounded on any time interval T. Recalling that the 
deterministic parQ of Up is a C 1,2l/ (R+, K d ) function, we can then use Lemma |2~51 
and Lemma l2.6l above, on the random part of Up for d = 1, 2, 3 to straightforwardly 
get the desired local Holder regularity for the direct solution of e ^ I f SLTBM ( a > u o), C^3 f 
as follows: we let q n = n + d for n G {0, 1, . . .} and let n = m + d for m = {0, 1, . . .}, 
we then have from Lemma 12.51 and Lemma 12.61 that 



E\Up(t,x) - Up{t,y)\ 2n+2d < C d \x- yf n+2d ^ , 

(3.7) < (2v-d){m+2d) 

E \Up(t, x) - Up(r, x)\ 2m+M <C\t- r\ ~» . 

for d = 1, 2, 3. Thus as in Theorem 2.8 p. 53 and Problem 2.9 p. 55 in [39] we get 
that the spatial Holder exponent is 7 S 6 ^0, -^-^pfg— ) and the temporal expo- 
nent is 74 g ^0, m ( 1 ~ rf /^H4 1 - d / t/ ) j y TO) n Taking the limits as to, n — > oo, we get 
7t € (0, and 7 S e (0, a d ), for d = 1, 2, 3. The proof is complete. □ 



4. Proof of the second main Theorem 

4.1. Regularity and tightness without the Lipschitz condition. As we men- 
tioned in Section T2.21 the finiteness assumption of Mp tq (t) and Mp tq (t) on T in 
Lemma 12.51 and Lemma 12.61 is for convenience only. We now proceed to show how 
to remove that assumption by showing it automatically holds under the weaker 
conditions QNLip| ). It is easily seen that if a is bounded then, for all spatial di- 
mensions d = 1, 2, 3, Mp tq is bounded on any compact time interval T = [0, T] (sec 
Remark 14.11 below) . The following Proposition gives an exponential upper bound 
on the growth of Mp_ q in time in all d — 1,2,3 under the conditions in ( |NLip[ ). 
The same result holds for Mp iq with only notational and obvious changes to the 
following proofs. 

Proposition 4.1 (Exponential bound for Mp, q ). Assume that Ug(t) is a solution of 
the B-ISLTRW SIE e SIE (a, un.n) onTx X d , and assume that the conditions in 

\_ l3-ISLTRW y ' u ' ' ' 

( |NLip[ ) are satisfied. There exists a constant C depending only on q, m&x x \uq(x)\, 
the dimension d, f3, and T such that 

M P ,S) < C (l + M p , q (a)ckJ ; < t < T, 

for every q> l,/3 <E {^;k € N} and d = 1,2,3. Hence, Mp, q {t) < Cexp{Ci} for < 
t < T,q > l,f3 G {^r;fcG N} , and d = 1,2, 3. In particular, Mp >q is bounded on T 
for all q> 1, p £ {^;fc e N}, andd= 1,2,3. 

The proof of Proposition |4~T1 proceeds via the following lemma and its corollary. 



Of course, the deterministic part of e^ IE (a, un) is, as discussed before, the integral 

fat* K Tx,y Afl u o(y) d V< and the random part is f Rd JjK™.^ a(U p(s,y))W (ds X dy). 
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Lemma 4.1. Under the same assumptions as in Proposition \4-l\ there exists a 
constant C depending only on q, max x |mq(x)| ; the dimension d, ft, and T such that 



Mp, q (t) < 




< t < T, 
t = 0, 

for every q > 1, f3 S k € N}, and d = 1, 2, 3. 



Proof. Fix (7 > 1, let Up D {t) = ^ K^^^Uo^) (the deterministic part of Up). 

Then, for any (t, i)£Tx X d , we apply Burkholder inequality to the random term 
&l R {t) to get 

2„ 



E 



2<l 



E 



(4.1) 



< C E 



E 



a<M 2a 2 (t/!( S )) 



5 d 



ds 



2<l 



Now, for a fixed point (t, x) € T x X d let /if be the measure on [0, t] x X rf defined 
(kS.^) /S d ds, and let |^| = Mt([M x Then, we can 



by dfj%{s,y) = 
rewrite (14.11) as 



(4.2) E U%{t) 



2q 



< c 



[0,t]xX d 



a 2 (C/|( S )) 



d^t(s,y) 



\tf\ q + \Uln(t)\ 2q 



Observing that fif/\[j^\ is a probability measure, we apply Jensen's inequality, the 
growth condition on a in QNLipD , and other elementary inequalities to (|4.2j) to 
obtain 



E 



2<l 



< C E 



< C 



= c 



[0,t]xX d 

E 



1 +E 



[0,t]xX d 



t—s:x,y 



a(E/*(a)) 



2<? dfi?(s,y) 
Im?I 



2g 



2<A 






2q 


J dtf{s,y) 


\ti\ q ~ l + c 


U% D {t) 





1 + E 



2</ 



rf.S 



I/'-/. 



i?- 1 



2g 



Using Lemma [12] we see that \fj,f\ is uniformly bounded for t < T and d = 1, 2, 3. 
So, using the boundedness of uq, and hence of U% D (t) by the simple fact that 
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J2 y ex d ^tvc% Ap = •*-! Lemma |2"T21 and the definition of Mp A (s), we get 




Here, R\ holds for d — 1,2,3. This implies that 

m,,S)<c(i + f^rds). 

Of course, Ma i9 (0) = sup^, |wo(a;)| 2 ' 3 < C, by the boundedness and nonrandomness 
assumptions on uq(x) in ( |NLip[ ). The proof is complete. □ 



Remark 4.1. It is clear that for a bounded a, Mp q is locally bounded in time. 
This follows immediately from Lemma [2~2l along with (|4.2[) above. 

Corollary 4.1. Under the same assumptions as those in Proposition \4-l\ there 
exists a constant C depending only on q, maxj |zt (a;) | , the dimension d, ft, and T 
such that 

M 0<q (t) < C ^1 + M 0>g (s)ds\ ,0 < t < T,q > l,P e {^;k G N} and d = 1,2,3; 
and hence 

Mp, q (t) < Cexp{C0; V0 < t < T, q > 1,0 G {^; k G N} , and d= 1,2,3. 



Proof. Iterating the bound in Lemma 14.11 once, and changing the order of inte- 
gration, we obtain 



(4.3) 



< C { 1 + C 



ds 



< G I 1 + / Mp, q {s)ds 



- + / Mp, q {r) 



ds 



(t - s)2i/ ( s - r)2v 



dr 



for d — 1,2,3. The proof of the last statement is a straightforward application of 
Gronwall's lemma to (14.31) . This finishes the proof of Corollary 14.11 and thus of 
Proposition ^. II □ 



The regularity, tightness, and weak limit conclusions for the /3-ISLTRW SIEs 
now follow. 
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Lemma 4.2 (Regularity and tightness). Assume that the conditions ( |NLip[ ) hold, 
and that < Ug (t) > is a sequence of spatially-linearly-interpolated solutions to 

the /3-ISLTRW SIEs h s J®(a,u ,n)} in (fi~23)) . ITien 

(a) -For every n, Up n (t) is continuous on M+ x R d . Moreover, with probability 
one, the continuous map (t,x) i— > Ug n (t) is locally ^t-Holder continuous in 
time with 74 e (0, ^=4) /or d = 1, 2, 3. 

(b) There is a /3-ISLTRW SIE weak limit solution to e g I ® SLTBM ( a i u o)} caZZ it 
Up, such that Up(t,x) is L p (f2,P) -bounded on T x R d /or every p > 2 and 

( 2v-d \ ~ - 

G H V 4,y / (T x M d ; R) for every d = 1,2,3 and ad G Id, w/iere ay 
and Id are as in Lemma \2.4\ 

Remark 4.2. Of course in part (a) above, even without linear interpolation in 
space, Ug(t) is locally Holder continuous in time with Holder exponent jt G (0, 2l ^ d ) 
ford =1,2, 3. 



Proof. For each n, let U% n (t) = Ug n D (t) + Ug n R (t) be the decomposition of 
Up n (t) in (|1.23[) into its deterministic and random parts, respectively. 

(a) By Lemma 11.21 Up n D (t) is clearly smooth in time; so it is enough to 
consider the random term Ug n R {t). We let q m = m + 2 for to e {0,1,.. .}, 



we then have from Lemma LL6I that 



(4-4) E ThM) - mJr) 



4+2m (2v-d)(m+2) 

<C\t-r\ z» 



for d — 1,2,3. Thus as in Theorem 2.8 p. 53 [39] we get that j t g 
(0, m ( 1 ~ d / 2v )+2- d - d /v j f or every to. Taking the limit as m oo, we 



2m+4 

get 7* G (0,^) ford = 1,2, 3. 
(b) By Lemma 12.11 it follows that Up n D (t) converges pointwise to the deter- 
ministic part of e^_f SLTBM (a,n ) in (JO); i.e., 

(4-5) lim Ul ntD (t) = [ K B t ^ y A "u ( y )dy. 



We also conclude from Lemma 12.51 and Lemma 12.61 that the sequence 
{U% R {t)\ is tight on C(T x R d ) for d = 1,2,3. Thus there ex- 

ists a weakly convergent subsequence I Up nk \ and hence a /3-ISLTRW 

I ' J k£N 

SIE weak limit solution U to e SIE (a.un). Then, following Skorokhod, 

/3-ISLTBM v ' u/ ' ° ' 

we construct processed Yp,k = Up,n k on some filtered probability space 
(£l s ,,^ s , {^f},P s ) such that with probability 1, as k -> oo, Yp, k (t, x) con- 
verges to a random field 1^(£, x) uniformly on compact subsets of T x R d 
for d = 1,2,3. Now, for the /3TSLTRW SIEs limit regularity assertions, 
clearly the deterministic term on the right hand side of (|4.5j) is C 1,2 " and 



2 ^As usual, = denotes equal in law or distribution. 
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bounded as in [3], so we use Proposition l4.il Lemma [2T5l and LemmaU 
to obtain the regularity results for the random part. We provide the steps 

here for completeness. First, Ypjt — Up. nk and so Proposition 14. II gives us, 
for each p > 2: 



(4.6) ElYp^x)] 1 



E 



U$,n k (*) < C < oo; V(t, x, k) g T x R d x N, d = 1, 2, 3, 



for some constant C that is independent of k,t,x but that depends on 
the dimension d. It follows that, for each (t,x) E T x R d the sequence 
{|Yfe(t, a;)| p } fe is uniformly integrable for each p > 2 and each d = 1,2,3. 
Thus, 

(4.7) E\Ufs{t,x)\ p = E\Yp(t,x)\ p = lim E \Yp tk (t, x)\ p < C < oo;V(t,x) € TxM d , 



(4.8) 



for all d = 1,2,3 and p > 2. Equation (|4.7|) establishes the L p boundedness 
assertion. In addition, for q > 1 and d = 1, 2, 3 we have by Proposition 14. II 

E \Y Pik (t, x)\ 2q + E \Yp 3k (t, y)\ 2q + E |K fljfc (r, x)| 2? 



< C 



< C; V(fc, r, i,i,i/)eNxT 2 x M 2 . 

So, for each (r, t, x, y) G T 2 xR 2 . the sequences jlYjg^t, x) — i^,fe(f, y)| 2? | 

and ||ys,fc(i, x) — Yp^(r, x)| 2? | are uniformly integrable, for each q > 1. 
Therefore, using Lemma 12.51 and Lemma 12.61 we obtain 

'E\Up(t,x) - U f} (t,y)\ 2q = E\Yp(t,x) - Y p (t,y)\ 2q 
= lim E \Y p , k (t,x) - Y 0tk (t, y)\ 2q < C d \x - y| 2 ^; a d g J d) 

,4 ' !)i < E|t^(t.«) ~ Up(r,x)\ 2q =E\Y (t,x) -Y (r,x)\ 2q 

= lim E\Y , k {t,x)-Y 0tk (r,x)\ 2q <C\t-r\ \ 

k — ^oo 

for d = 1,2,..., 3. The local Holder regularity is then obtained using 
exactly the same steps as in Q3.7P and the following conclusions. 

The proof is complete □ 



4.2. Recalling the K-martingale approach. For the article to be self-contained, 
we now recall and briefly discuss the K-martingale approach from [1] — adapting 
it to this paper's settinepl This approach is tailor-made for kernel SIEs like 
e^ IE „ (a.un) and other mild formulations for many SPDEs on the lattice. The 

/3-ISLTBM v ' u/ J 

first step is to truncate to a finite lattice model as in (II -25[) . Of course, even after we 
truncate the lattice, a remaining hurdle to applying a martingale problem approach 
is that the finite sum of stochastic integrals in fjl .25[> is not a local martingale. So, 
we introduce a key ingredient in this K-martingale method: the auxiliary problem 
associated with the truncated /3-ISLTRW SIE in ()1.25p . which we now give. Fix 



30 A11 we need to adapt it here is a notational change, replacing the BTRW transition density 
in [1] with the /3-ISLTRW one. 
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(l,n) G N 2 and r G M.+ . We define the r-auxiliary ,3-ISLTRW SIE associated with 

(|L25l) on [0, t] x X* by 

(Aux) 



x€X d n \ X* 



> o, 



where the independent BMs sequence {W%} yGX d in (IAux[) is the same for all r 
as well asieX^j. We denote (jAuxft by e ^g^ m ( a > u o, ™, t). We say that the 

pair of families ({^ln,l} T>Q » TC} v6 x* ,) solves {e^ E RW (a, u , n, I, r)}^ on 

a filtered probability space (fi, J?", {^}, P) if there is one family of independent BMs 
(up to indistinguishability) {W»(t); < i < oo} ygX<1 ; on (0, J?, { J^}, P) such that, 
for every fixed r G M+ 

(a) the process | i(i),^(;0<f<T,i£ XjJ j has continuous sample paths 

in f for each fixed a; G X^ and Xl'^ t (t) G for all x G X^ for every 
< i < r; and 

(b) equation (jAuxft holds on [0,r] x X^, P-almost surely. 

Naturally, implicit in our definition above the assumption that, for each fixed r G 



we have 



= 1; Vas,yeX^,,0<t<T. 



For simplicity, we will sometimes say that X^ n l = {aJ* ; (i),^;0 < t < t,x G 

is a solution to (|Aux[) to mean the above. Clearly, if Al'^ satisfies (IAux[) then 

Up t n t i{T) := X T p X n ; (r) satisfies (|1 .25[) at i = r for all a; G X„. Also, for each n and 
each d = 1, 2, 3 



S"— s;x,y 



5>, 



< 



cd/2 



In addition, for each fixed r G ffi+ and each fixed x,y £ X^ , we have for a solution 
A^ n / to (|Aux)l that 



G & s 



Vs < r, 



since, of course the deterministic ¥^_ s s n J. A y / 'S?/ 2 G and a(XZ' v n ; (s)) G ,^ s . Thus, 
if XJ] n l solves (|Aux[) ; then, for each fixed r > and x, y G X^ ; , each stochastic 
integral in (|Aux[) 



is a continuous local martingale in t on [0, r]. This is clear since by a standard 
localization argument we may assume the boundedness of a (|a(u)| < C); in this 
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case we have for each fixed x, y £ X„ ; and r G R+ that 



E 



T T ' x - y (t) 



- x Ls,r ^oo = -EOT, 



r < t < T. 



So, the finite sum over X^ l in (|Aux|) is also a continuous local martingale in t on 
[0, t]. I.e., for each r > and iEXJ, 



M l%,m--= E r«;;L, T tew)^(^^;o<t<rU 



//c.loc 
*2 



with quadratic variation 



(is 



where we have used the independence of the BMs {W^} weX d within the lattice X^ ; . 

For each r > 0, we call M x ^ r n l a kernel local martingale (or K-local martingale). 

There is another complicating factor in formulating our K-martingale problem 
approach that is not present in the standard SDEs setting. To easily extract so- 
lutions to the truncated /3-ISLTRW SIEs in (|1.25j) from the family of auxiliary 

problems | e ^ u f S L^ l w'' a ' w °' n ' r )} m ^ UX P ■ we wan t the independent BMs se- 
quence {W^} ygX d to not depend on the choices of t and x. I.e., we want all the 

K-local martingales in (|Aux[) to be stochastic integrals with respect to the same 
sequence {W y } yeX d , regardless of r and x. With this in mind, we now formulate 

the K-martingale problem associated with the auxiliary /3-ISLTRW SIEs in (jAuxl) . 
Let 

(4.11) C n .i := {it : M+ x (Xf M ) 2 ->■ R 2 ;t \-t u Xl > X2 (t) is continuous Vx 1: x 2 \ ■ 

For u € C n j, let u Xl ' X2 (t) = (u Xl (t) , u x 2 2 (t)) with u x (t) = u x ' x (t); and for any 
T\ , t~2 > and any xi,X2,y S X^ l let 



(4.12) 



,d/2 



^a(u?(t))- 



-t;xj,y 



,d/2 



a{u y At))- 1<*,J<2, 



(we are allowing the cases ri = T2 and/or xi = #2) where for typesetting con- 
venience we denoted the points (t^t^) and (xj,Xj) by r, ( j and a^j, respectively. 
We denote by <9i and 9?- the first order partial derivative with respect to the i-th 
argument and the second order partials with respect to the i and j arguments, 
respectively. Let C 2 = C 2 (M 2 ;!L) be the class of twice continuously differentiable 
real-valued functions on R 2 and let 



(4.13) C 2 = {/ <E C 2 ; / and its derivatives up to second order are bounded} 
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Now, for ti,t 2 > 0, for / G C§, and for (t,xi,x 2 ,u) G [0, n A r 2 ] x (x^^ x C n> / 
let 



l<i<2 

+ ^ E ^/(« Bi,Bj (*)) E T ^ («"(*)) 

l<i,j<2 



(4.14) 



Let Xg n j = j-X^* i(t);0 < t < r, x G be a continuous in i adapted real-valued 
process on a filtered probability space (Q, J^", {^i},P). For every ti,t 2 > define 
the two-dimensional stochastic process ZV"*^ 

(4-15) {z;%T 2 (t) = (x™(t),x%%(t)) ; (t,*i,z 2 ) e [o,n a r 2 ] x (x^) 2 } 

with Z*;^ a (t) = (X^(t),^(t)) and let f/*^ = M*i), u (x 2 )). We say 
that the family i XI n l > satisfies the K-martingale problem associated with the 

auxiliary /3-ISLTRW SIEs in (Tftux|) on R+ x X^ if for every / G Cg, < n, r 2 < oo, 

t = n A t 2 , f G [0, r], xi, x 2 G Xfj j, and x G XjJ \ X^ ; we have 
(KM) 

/(^T 2 (*)) - /(^o IlX2 ) - J* K 1 ' 3 /) * a , G ^ oc ; 

^(*) = ^,».i>(*)- 
We are now ready to state the equivalence of the K-martingale problem in fKMj 
to the auxiliary SIEs in (jAuxj) and its implication for the /3-ISLTRW SIE in JQ5]) . 
This result is of independent interest and is stated as the following theorem 31 !. 

Theorem 4.1. The existence of a solution pair ^\x^ nl ^ , {W^} yeXd \ to 

I e aux-s'iE ' ( a , uq, n, I, t) \ in (jAuxj) on a filtered probability space (O, JP, {^ t }, P) 

^ p lOLilJXW J ^_|v,Q 

is equivalent to the existence of a family of processes < XZ n l > satisfying (|KMj) . 

I ^' ' J r>0 

Furthermore, if there is 

{ X kn.i} >0 satisfying (|KM]l th en there is a solution to 

(fl~25l) on R+ x X*. 

The proof follows the exact same steps as the proof of Theorem 1.3 in [TJ and 
will not be repeated. 



4.3. Completing the proof of the second main result. We now complete the 
proof of Theorem 1 1.31 In Section [2?2l and Section |4~T1 we assumed the existence of a 
/3-ISLTRW SIE solution and we obtained regularity and tightness for the sequence 
of lattice SIEs < e s Jf, r ^, T (a,Uo,n) } . This, in turn, implied the existence and 

I P- lbLiRW J riGN* 

regularity for a /3-ISLTRW SIE limit solution to our e ^f SLTBM ( a , w o) in (|1.14|) . To 



31 This is because it is easily adaptable to many mild formulations of SPDEs, of different orders, 
not just for the /3-ISLTBM SIEs. Since we don't prove uniqueness under less than Lipschitz 
conditions for our /3-ISLTBM SIEs, we have not explicitly mentioned the uniqueness implications 
of our K-martingale approach. More on that in future articles. 



/3-INVER.SE-STABLE-LEVY-TIME BROWNIAN MOTION SIES ON K + x R d 



35 



complete the existence of the desired double limit solutiorrl for e SIE (a.uo) 

^ — /3-ISLTBM v ' u; 

it suffices then to prove the existence of a solution to e SIE (a.un.n) for each 

^ /3-ISLTRW v ' u ' ' 

fixed n £ N*, under the condition ( |NLip[ ), that is uniformly L p (Sl,F) bounded on 
[0, T] x X d for every T > and every p > 2. We establish this existence via the 
K-martingale approach just recalled and adapted from pQ, using Theorem 14. II 

First, the following proposition summarizes the results in this case for the (3- 
ISLTRW SIEs spatial lattice scale0. 

Proposition 4.2 (Existence for /3-ISLTRW SIEs with non-Lipschitz a). Assume 
the conditions ( |NLip[ ) hold. Then, 

(a) For every (n, I) £ N*xN, every /3 = 1/v £ {l/2 fe , k £ N}, and for every p > 
2, there exists an LP -bounded solution Ug n to the truncated /3-ISLTRW 
SIE (ll.25[) on TxX^. Moreover, if we linearly interpolate Ug n ; (t) in space; 
then, with probability one, the continuous map (t,x) t— > Ug nl (t) is locally 
It-Holder continuous in time with j t £ (0, ) for v — /3 _1 £ {2 k ; k £ N} 
and d = 1,2, 3. 

(b) For any fixed n £ N* , £/ie sequence | L^f n ;(t) | o/ linearly-interpolated 

solutions in (a) has a subsequential weak limit Up :1l in C(T x M d ;M). We 
thus have a limit solution Ua, n to eS g E ISLTBW { a i u 0i n )> an d Up tn is locally 74- 
Holder continuous in time with -y t € (0, for v = j3~ x £ {2 k ; k £ N} 

and d = 1,2, 3. 



Proof. 

(a) First, recall that the deterministic term U% D (t) in (|1.25l) is completely de- 
termined by uq- Moreover, under the conditions in ( |NLip[ ) on uq, f7| D (t) 
is clearly bounded and it is smooth in time as in Remark 1 1.1 1 Fix an arbi- 
trary T > 0, and let T = [0, T]. We now prove the existence of a family of 

adapted processes \ XZ , > satisfying our K-martingale problem (|KMI) , 
I ^' ' J i-eT 

which by Theorem 14.11 implies the existence of a solution to the Z-truncatcd 
/3-ISLTRW SIE (Ti~25l) on T x x£. On a probability space (St, =T, {^},P) 
we prepare a family of r-independent BMs {W%(t)} y£X d ; . For each r e T 

and each i = 1,2,... define a continuous process Xl nli on [0,r] x X„ 
inductively for k/2 l < t < ((k + l)/2*) At (k = 0,1,2,...) as follows: 
X p%,l,i(°) = u o( x ) G X «) and if X^ l t (t) is defined for t < k/2\ then 



^ 2 The type of our lattice limit solution to e^ IE (a, no) m 11.141 depends on the conditions: 

/3-IoLTBM 

under the Lipschitz conditions dLipfc we get a direct solution to the lattice SIE e| IE „„,,, (a,uo,n) 

' ' /3-ISLTRW 

for every n and a direct /3-ISLTRW SIE limit solution to e^ IE „ (a, tin) (see Theorem IAJJ ; 

J ' /3-ISLTBM v ' ' v 1 " 

whereas under the non-Lipschitz conditions in | |NLip[ | we obtain a limit /3-ISLTRW SIE so- 
lution, thanks to our K-martingale approach, and a /3-ISLTRW SIEs double limit solution to 

^-fsLTBM^'" )- 

■^We remind the reader that we will, without further notice, suppress the dependence on 
whenever it is more convenient notationally to do so. 
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we 



define Xf* nU {t) for k/2 l <t<{{k + l)/2 4 ) A r, by 



(4.16) 



Y X ' T ( k 



Ul n , D {t) - U* D (£) 



X G 



where A j, W» = - Clearly, X^„ M is the solution to the 



equation 



(4-17) j E r^ W ,r(^ ,T (^(*)))^( a ) + ^ ( n,Ij(*); ^ eX lf 



with X^ T (0) = u (a:), where = fc/2 l for fc/2 4 < * < (fc + 1)/2 J A 

(fc = 0,l,2,...). 

2q 



Now, for 9 > 1, let M£ , ,(<) = sup xeX , E X^ l t {t) 



By the bound- 



(4.18) 



edness of U% n D (t) over the whole infinite lattice X^ , we have 

2q 



Mp*,i,i{t)<C+ sup E 



X 



(*) 



(4.19) 



Then, replacing X^ by X^ ; and following the same steps as in the proof of 
Proposition 14. 1[ we get that 

sup sup Ml ,At)<C, d=l,2,3, 

rSTte[0,r] 

where, here and in the remainder of the proof, the constant C depends only 
on q, (3, maxj |uo(x)|, the spatial dimension d — 1,2,3, and T but may 
change its value from one line to the next. Remembering that 5 n \, as 
n /*• oo and ii e N*, the independence in I is trivially seen since Lemma 
implies 



E 



H;x,y 



< 



E 



t:x,y 



6 n 



< 



c 

ld/2h 



:Vd = 1,2,3, Z e N 



Similarly, letting X^ l 'J i R denote the random part of X x ^ T n l i on the trun- 



cated lattice X^ t , using (|4. 19[) . and repeating the arguments in Lemma l2~5l 
and Lemma 12.61 — replacing X^ by j and noting that Lemma 12.31 and 
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Lemma 12.41 hold on X^ z — we obtain 



E 



2<1 



~ X | n^i flft) + E X fi',n,l,i,R^ ~ X f^,l,i,R^ 



(4.20) 



<C d \x-y\^ a ";a d Gl d , 

X P,n,l,i,R$) ~ X P,n,l,i,R( T ) 
(2v-d)q 

<C\t-r\ 2„ , 



2-/ 



E 



X p,n,l,i,R^) X p,n,li,R( r ) 



2q 



2q 



for all x,y G X^, r,t G [0,n A t 2 ], ri,r 2 G T, and d = 1,2,3. It fol- 



lows that, for every point tj.,2 = (71,72) G T 2 , there is a subsequence 
{{ X }]n,l,i m > X 2 ,n,l,i m )} _ t on a probability space (0 TM) # riiJ) P Tl|! ) such 



that (X£ n j. m ,X% n:l . 



uniformly on compact subsets of [0, t\ A r 2 ] x X^, asm->oo a.s. Let Tq = 
T n Q, where Q is the set of rationals, and define the product probability 
space 



vri, 2 eTg T!, 2 eTg n, 2 eK 



If s < t, then for every / G C 2 (R 2 ;R), n,r 2 G T Q \ {0}, t G [0,n A r 2 ], 
xi,a;2 G X^ j , and for every bounded continuous F : C (R + ;R 2 ) — > R that 
is £3 S (C (R+; R 2 )) := a (z(r); < r < s)-measurable function, we have 



E 5 



(4.21) 



{/TO n,a (*))-/TO n,a W) 

«-/) (r, xx.xa, ^)dr j F fef^G) 



= lim 1 

m— foo 
ft 



{/(^^"(*))-/(%^w) 



r,xi,a;2, Z 



where, by a standard localization argument, we have assumed that a is also 
bounded; and where Z^ 1 ^ z and Zfp~£ l im are obtained from the definition 

Tl.2 



of Z;% t in 613 by replacing J^ n>l by X^ n l and X|„, Mm , j = 1,2, 
respectively. The operator s^^f is obtained from rfj 1,2 by replacing 

(^(0) in b y T S m (t).T 4lJ («"(&»(*)))■ Also, obviously, 

for any r G Tq and t G [0, r] 

(4.22) X; : T n l (t) = lim X|;; Mm (t) = tft B , D (i); x G X^ \ X? M , a.s. P. 
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(4.23) 



(4.24) 



(4.25) 



It follows from (|4.21l) and (|4.22[) that n ; | satisfies the K-martingale 
problem (|KM[) with respect to the filtration {J^*}, with 



fl a {^|,n,l(«)5 « < (^ + e) A r, r e T Q n (t, T] } 

e>0 



Thus, by Theorem 14. 1[ with t G R+ replaced by r G Tq, there is a solu- 
tion U x p n l {t) to the Z-truncated /3-ISLTRW SIE CTJ5]) on T Q x x «- Use 
continuous extension in time of C/f ra z (t) to extend its definition to T x X^, 
and denote the extension also by Ug nl (t). Clearly Up nl {t) solves the l- 
truncated /3-ISLTRW SIE (H~25l) on T x X^. 

2q 



Now, for q > 1, let Ms, q ,i(t) — sup xe; 
boundedness of Up nD (t), implies 



d E 



M M {t)<C+ sup E Ul n j(t) 



2<l 



As above, the 



Then, replacing X^ by X^ , and following the same steps as in the proof of 
Proposition 231 we get that 

Mp, q ,i(t)<C, Vi G T,/3 G {l/2 k ;k G N} and d = 1,2, 3. 

Similarly, letting t/| n ; denote the random part of Up n on the 
truncated lattice X^ t , using (|4.24l) , and repeating the arguments in Lemma| 
and Lemma 12.61 — replacing X^ by X^ ; and noting that the inequalities in 
Lemma 12.31 and Lemma 12.41 trivially hold if we replace X^ by X^ l — we 
obtain 



2<I 



<C d \x- 



E 



2<l 



<C\t 



y\ 2qad ; a d el d , 

(2v-d)q 



for all x,y € Xf u , r,t G T, and d = 1,2,3. By Remark [TTJ Up nD {t) is 
differentiable in t. So, linearly interpolating Up nl (t) in space and using 
(|4.25p and arguing as in the proof of part (a) of Lemma l4~2l we get that the 
continuous map (t, x) i— > Up nl (t) is locally 7t-H61der continuous in time 
with 7f G (0,^r) for v = yS - ' 1 G {2 fc ;A: G N} and d = 1,2,3. 



(b) Clearly, Ug D (t) in (jl.25l) is the same for every I, so it is enough to show 
convergence of the random part C/l t R {t). Using ()4.25|) we get tightness 
for {l7, 



n i flWj anc ^ consequently a subsequential weak limit l^g.n, which 
is our limit solution for e ^f SUTRW i a t uq, n). For the regularity assertion, 



p n D (t) is smooth and bounded as noted above. So, using (|4.24l) and 



(|4.25p , and imitating the argument in the proof of part (b) of Lemma 14.21 
(remembering that here we are taking the limit as I — > oo); we get the 
desired LP boundedness for Ug, n as m Proposition 14. II and the spatial and 
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(4.26) 



temporal moments bounds in Lemma 12.51 and Lemma 12.6 



2-/ 



< c 



2q 



Ut n , R (t) ~ Ul R (t) < C d \x - y| 2 ^; a d G I d , 

2q (2v-d)q 

Ul n , R (t)-Ul R (r) <C\t-r\ 2, , 



for (t, x, n) G T x X d x N* and for v — /3 _1 G {2 fc ; fc G N}, d = 1, 2, 3, and 
q > 1 and the desired Holder regularity follows. 

The proof is complete. □ 

We now get Theorem 11.31 for e|™ BM (a, uq) as the following corollary. 
Corollary 4.2. Theorem \l.&\ holds. 



Proof. The desired conclusion follows upon using the argument in the proof of 
part (b) of Lemma 14.21 along with Definition 11.41 and the L p -boundedness and the 

spatial and temporal moments bounds for < tJp, n \ that we got in (|4.26[) above. □ 



Appendix A. Limit solutions in the Lipschitz case 

We now state prove our lattice-limit solution existence, uniqueness, and regular- 
ity for our BTBM SIE on R+ x R d under Lipschitz conditions. 

Theorem A.l (Lattice- limits solutions: the Lipschitz case). Under the Lipschitz 
conditions there exists a unique- in-law direct /3-ISLTRW SIE weak-limit solution to 
e s B ™ BM (a,u ), U, such that U(t,x) is Lf(fi,P) -bounded on T x R d for every p>2 
(*=& A iV 

and U G HV iv ) V 2 ) (T x for every d — 1,2,3. 

Theorem IA.1I follows as a corollary to the results of Section 12.21 combined with 
the following proposition. 

Proposition A.l. Under the Lipschitz conditions ( |Lip[ ) there exists a unique direct 
solution to e ^ J f SLTflw X a : u 0i n ), Up^ n , on some filtered probability space (O, {^t}, P) 
that is LP (fi,P) -bounded on [0,T] x X d for every T > ; p > 2, n G N* ,and 
d = l,2,3. 

The proof of Proposition IA.1I follows the same steps as the non-discretization 
Picard-type direct proof of the corresponding part in the continuous case in Sec- 
tion [3l with obvious changes, and we leave the details to the interested reader. 

Corollary A.l. Theorem ] A. 1\ holds. 



Proof. The conclusion follows from Proposition lA.il Lemma |2~51 Lemma |2~51 and 
Lemma l4~2l (b). □ 
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Remark A.l. With extra work, it is possible to prove the existence of a strong 
limit solution under Lipschitz conditions. We plan to address that in a future 
article. 

Appendix B. Glossary of frequently used acronyms and notations 
I. Acronyms 

• BM: Brownian motion 

• BTBM: Brownian-time Brownian motion. 

• BTBM SIE: Brownian-time Brownian motion stochastic integral equa- 
tion. 

• BTP: Brownian-time processe. 

• BTP SIE: Brownian-time process stochastic integral equation. 

• BTC: Brownian-time chain. 

• BTRW: Brownian-time random walk. 

• /3TSLTRW DDE: Brownian-time random walk differential-difference 
equation. 

• /3TSLTRW SIE: Brownian-time random walk stochastic integral equa- 
tion. 

• DDE: Differential difference equation. 

• KS: Kuramoto-Sivashinsky. 

• RW: Random walk. 

• SIE: Stochastic integral equation. 

II. Notations 

• N: The usual set of natural numbers {1, 2,3,.. .}. 

• K t . x y : The d-dimensional continuous-time random walk transition 
density, starting at x g X^ and going to y £ X^ in time t. 

• Kftxy'- Tne density of a d-dimensional BM. 

• Kf.Q 1 /. The density of a 1-dimensional BM, starting at 0. 

• 1 : The kernel or density of a e?-dimensional Brownian-time Brow- 
nian motion. 

• K^. x S y' fi : The kernel or density of a c£-dimensional Brownian-time 
random walk on a spatial lattice with step size S n in each of the d- 
dimensions. 

• e SIE (a, m ): The BTBM SIE with diffusion coefficient a and initial 

BTBM v ' v/ 

function uq. 

• e s a l t Trj ,(a, uo,n): The /3-ISLTRW SIE on the lattice X^ = S n Z d with 
diffusion coefficient a and initial function Uq. 
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